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SECTION  I 
INTRODUCTION 


1.  OBJECTIVES 

In  designing  the  Air  Force  Refuse-Collection  Scheduling  Program  (RCSP), 
the  fundamental  objective  was  to  reduce  collection  costs.  The  most  signifi- 
cant cost  reduction  is  effected  by  a reduction  in  the  number  of  collection 
trips  used  to  service  a given  region.  If  a collection  crew  can  be  dropped 
from  the  fleet,  the  cost  of  manpower  will  be  cut.  In  addition,  fuel  and 
maintenance  costs  will  be  lessened  if  the  total  mileage  traveled  by  the 
collection  fleet  can  be  reduced.  The  first  objective,  then,  is  to  generate 
a collection  schedule  that  calls  for  the  theoretical  minimum  number  of  trips. 

The  second  objective  was  to  produce  a desirable  collection  schedule  that 
would  be  relatively  easy  to  implement.  Some  of  the  features  considered  de- 
sirable for  a refuse-collection  trip  by  one  vehicle  are  the  use  of  as  few  U- 
turns  as  possible,  collection  of  all  refuse  on  a block  by  one  pass  of  the 
vehicle,  spacial  clustering  of  all  streets  serviced  by  the  vehicle,  and  mini- 
mization of  the  time  and  distance  required. 

A third  objective  was  to  produce  a computer  program  that  is  easy  to  im- 
plement and  runs  quickly.  The  requirements  for  easy  implementation  include 
a simple  technique  for  describing  the  collection  region  to  the  computer,  no 
need  for  human  decisions  during  the  program  execution,  and  computer  output 
in  forms  that  can  be  used  immediately.  It  was  desirable  to  develop  the  pro- 
gram at  a cost  less  than  that  of  any  other  scheduling  method  and  in  a short 
time.  The  cost-performance  balance  finally  used  was  that  of  coming  as  close 
as  possible  to  the  performance  goals  using  roughly  a one-man-year  effort 
over  two  calendar  years. 

The  performance  finally  achieved  comes  very  close  to  the  stated  objectives. 
RCSP  usually  achieves  the  theoretical  minimum  number  of  trips  without  human 
intervention  and  always  achieves  the  minimum  with  a very  little  human  in- 
tervention. Implementation  is  fairly  straightforward,  although  the  desir- 
ability of  the  final  schedule  depends  on  the  care  with  which  the  map  of  the 


collection  region  is  reduced  to  computer  input  form.  Provisions  exist  for  hu- 
man improvement  of  the  computer-generated  routing  before  the  final  maps  and 
schedules  are  produced.  Spacial  clustering  of  the  streets  in  the  collection 
region  usually  is  very  good  for  all  but  one  trip. 

2.  SCOPE 

This  section  (Volume  I)  of  the  report  describes  the  first  program,  RCINPT. 
A program  overview  is  given,  followed  by  a thorough  description  of  the  logic 
involved  in  the  map  processing.  A skeleton  of  the  logic  flow  is  provided.  In- 
put and  output  files  are  described.  Program  requirements  and  restrictions, 
error  messages  and  error-handling  techniques,  definitions  of  important  symbols, 
and  a running  time  estimate  are  also  provided. 


SECTION  II 
PROGRAM  OVERVIEW 


Program  RCINPT  serves  two  purposes:  it  plots  the  map  input  data  to  verify 
their  accuracy,  and  it  determines  the  total  amount  of  refuse  for  which  collec- 
tion is  to  be  scheduled.  Since  the  refuse-quantity  computation  is  embedded  in 
the  map  processing  procedure,  the  functional  description  of  RCINPT  will  stress 
the  latter. 

RCINPT  receives  three  types  of  input.  The  first  data  record  contains 
street  names.  (The  program  provides  printouts  of  these  names  that  can  be 
checked  for  accuracy;  the  data  are  then  saved  on  disk  for  use  by  program 
PHASE4. ) The  second  data  record  describes  maps  to  be  plotted  by  RCINPT  as 
an  aid  to  debugging  the  input  description.  The  map-description  cards  are 
in  free  format,  with  refuse-quantity  information  embedded  in  the  map  de- 
scription. 

The  program  consists  of  a main  program  named  RCINPT  and  10  nonsystem 
subroutines.  RCINPT  calls  subroutine  STRINP,  which  reads  the  problem  title 
and  the  street-name  data.  STRINP  prints  the  title,  the  street  numbers,  and 
the  street  names.  The  street-name  information  is  buffered  out  to  file  TAPE3. 
Control  returns  to  RCINPT. 

RCINPT  reads  the  second  data  record,  which  consists  of  the  bounding  co- 
ordinates and  the  sizes  of  maps  to  be  plotted  by  subroutine  MAPPLT.  If  any 
maps  are  to  be  plotted,  RCINPT  calls  subroutine  MAPGRID,  which  draws  the  axes 
and  a grid  on  which  the  map  will  be  plotted.  MAPGRID  uses  subroutine  AXIS  to 
draw  each  axis  and  to  append  tic  marks  to  the  plot.  Subroutine  AXIS  uses  sub- 
routine NUMBER  to  append  numbers  to  the  axis. 

After  drawing  the  grid,  RCINPT  reads  two  cards  from  the  first  map  descrip- 
tion, which  give  various  parameters  pertaining  to  the  map.  These  parameters 
are  then  printed  out.  The  remainder  of  the  map-description  data  consist  of 
groups  of  strings  describing  street  connections  and  positions.  Refuse- 
quantity  information  is  embedded  in  these  strings.  The  strings  are  read  in  free 
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format,  using  subroutine  LINEIN.  The  strings  are  read  in  groups  corresponding 
to  street  segments;  as  each  group  of  numbers  is  read,  the  information  is  stored. 
Subroutine  M0VE5  and  function  IF IND  are  used  to  position  the  data  in  the  corre- 
sponding arrays.  Each  string  is  terminated  by  coordinates  and  shape  informa- 
tion. The  shape  information  is  processed  by  subroutine  SHAPCOM.  Coordinates 
of  points  on  each  street  segment  are  determined  by  subroutine  COORD. 

When  each  string  has  been  completely  read  in  and  processed,  subroutine 
MAPPLT  is  called  to  plot  the  information  from  that  string.  This  process  is  re- 
peated until  all  strings  in  the  record  have  been  processed.  Additional  records 
of  map  data  are  processed  in  the  same  manner,  and  the  reading  of  map  data  term- 
inates at  an  end-of-file  or  two  consecutive  end-of-record  cards. 

When  all  of  the  map-description  data  have  been  read,  the  information  about 
each  street  segment  is  printed  out.  Information  about  street  segments  is 
written  on  disk  file  TAPE1 , and  information  about  intersections  (nodes)  is 
written  on  disk  file  TAPE2.  If  more  than  one  output  map  is  requested,  RCINPT 
calls  subroutines  MAPGRID  and  MAPPLT  repeatedly  until  all  output  maps  have  been 
completed. 

The  flow  of  control  from  one  subprogram  to  another  is  shown  in  Figure  1. 
Within  each  subprogram,  only  the  first  call  to  each  other  subprogram  is  shown. 
(Four  of  the  subroutines  shown  in  Figure  1--PL0TS,  SYMBOL,  PLOT,  and  EXIT--are 
subroutines  from  the  basic  Calcomp  software  package  and  are  not  included  in 
the  description  of  program  RCINPT.) 
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SECTION  III 
PROGRAM  LOGIC 


The  logic  for  program  RCINPT  is  described  from  three  viewpoints.  The 
first  description  is  task  oriented.  The  second  view  is  data-storage  oriented 
and  includes  a discussion  of  the  preparation  of  data  for  use  by  subsequent 
programs.  The  third  view  describes  each  subroutine  in  terms  of  its  purpose 
and  the  manipulations  performed  within  it. 

In  the  descriptions  that  follow,  certain  terms  are  used  to  represent  por- 
tions of  arrays.  Assume  that  array  X has  permissible  subscripts  of  from  1 to 
100.  The  terms  low,  first,  and  front  refer  to  data  at  or  near  X(l).  The  terms 
high,  last,  and  end  refer  to  data  at  or  near  X ( 1 00 ) . 

1.  PROGRAM  TASKS 

Program  RCINPT  serves  two  primary  tasks:  it  reads  the  map-description 
data,  and  it  plots  the  output  maps. 

Program  execution  begins  in  main  program  RCINPT.  The  Calcomp  plot  pack- 
age is  initialized  by  a call  to  subroutine  PLOTS.  Program  RCINPT  calls  sub- 
routine STRINP,  which  reads  in  the  problem  title,  the  street  numbers,  and  the 
names  associated  with  those  numbers. 

The  main  program  then  reads  the  second  record  of  data,  which  will  be  the 
descriptions  of  the  regions  to  be  plotted  on  maps.  From  1 to  10  cards  must  be 
present  in  this  record.  The  description  of  the  collection  region  may  consist 
of  several  maps,  each  producing  a record  of  map-description  data.  The  first 
two  cards  of  each  map-description  record  consist  of  scale  factors  and  default 
parameters.  These  are  read  by  the  main  program  using  formatted  read  state- 
ments. The  problem  title  and  the  map  parameters  are  printed.  The  remaining 

data  cards  in  the  map-description  record  consist  of  free-format  strings  that 
describe  the  street  connections,  the  refuse  quantities,  and  the  street  geom- 
etries. These  cards  are  read  by  several  calls  to  subroutine  LINEIN. 


i 


The  first  call  occurs  at  statement  210  in  the  main  program.  The  street 
number  and  the  number  of  the  first  node  on  the  card  are  read  by  this  call  to 
LINEIN.  The  next  four  items  on  the  map-description  string  may  be  repeated  many 
times.  These  four  i tems--street-segment  length,  houses  on  right  side,  houses 
on  left  side,  and  end  node  number--are  read  <by  the  call  to  LINEIN  at  statement 
250  in  the  main  program.  The  data  obtained  on  this  call  are  processed  and 
stored  in  the  segment  data  table  and  the  node  data  table.  Control  then  returns 
to  the  call  to  LINEIN  at  statement  250  until  either  the  end  of  the  card  or  a 
slash  or  left  parenthesis  is  encountered.  If  a slash  is  encountered,  a speed 
limit,  a one-way  indicator,  the  number  of  sides  collected  on  one  pass,  and  a 
refuse  quantity  adjustment  factor  may  follow.  These  are  read  by  the  call  to 
LINEIN  preceding  statement  280  in  the  main  program. 

After  these  data  have  been  processed,  LINEIN  is  called  at  statement  310 
to  obtain  the  coordinates  of  the  first  node  and  the  shape  code.  Following 
statement  320,  LINEIN  is  called  to  obtain  the  coordinates  of  the  final  node 
of  the  string.  This  completes  the  reading  of  a map-description  string.  The 
coordinates  of  the  initial  and  final  nodes  of  the  string  are  processed  accord- 
ing to  the  mode  of  coordinate  use  specified  on  the  first  data  card  of  the  map- 
description  record.  If  average  coordinates  are  requested,  the  coordinates  are 
accumulated  into  a running  average  for  each  node.  If  the  last  coordinates  are 
to  be  used,  all  previous  values  are  overwritten.  If  only  the  first  specifica- 
tions of  coordinates  are  to  be  used,  no  subsequent  values  are  retained.  If  a 
new  map  with  a new  coordinate  system  in  a new  map-description  data  record  is 
being  processed,  the  translations  for  that  map  relative  to  the  initial  map 
description  are  found  by  using  either  the  initial  or  the  final  node  of  the 
first  map-description  string.  Therefore,  one  of  these  nodes  must  have  been 
specified  on  a previous  map. 

Processing  in  the  main  program  continues  by  the  assigning  of  an  individ- 
ual shape  code  to  each  segment  in  the  string.  Where  possible,  simplified  shape 
codes  are  assigned  to  the  segments.  For  example,  if  segments  on  a string  with 
an  angle  shape  code  exclude  the  vertex,  these  segments  are  assigned  a shape 
code  indicating  a straight  segment.  This  processing  takes  place  between  state- 
ments 600  and  690  in  the  main  program.  Where  collection  is  from  one  side  of 
the  street  at  a time,  the  segment  is  processed  twice,  one  way  in  each  direction. 
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Subroutine  MAPPLT  is  called  to  draw  the  contribution  to  the  first  output 
map  from  each  segment  in  the  string.  This  process  of  reading  strings  and 
plotting  their  contributions  to  the  first  map  is  repeated  until  all  of  the  in- 
put map-description  records  have  been  read.  When  the  first  output  map  has  been 
completed,  the  node  and  segment  data  are  written  to  disk  or  tape.  Subroutines 
MAPGRID  and  MAPPLT  are  called  repeatedly  to  draw  subsequent  maps.  Processing 
is  then  terminated  and  control  returns  to  the  system. 

As  each  string  in  the  first  map-description  record  is  processed,  subrou- 
tine MAPPLT  is  directed  to  draw  all  of  the  segments  in  that  string.  It  is  di- 
rected to  draw  all  of  the  segments  on  subsequent  calls  for  subsequent  maps. 

In  all  cases,  parts  of  the  map  outside  the  bounds  specified  on  the  cards  in 
the  second  data  record  are  not  drawn. 

Before  each  segment  is  drawn,  subroutine  SHAPCOM  is  called  to  set  up  the 
parameters  used  to  determine  the  position  on  the  segment  of  points  at  a given 
distance  from  the  start  of  that  segment.  The  actual  coordinates  of  points  on 
the  segment  are  returned  by  subroutine  COORD.  The  loop  from  statements  130  to 
170  in  subroutine  MAPPLT  generates  points  on  each  segment  and  plots  a line 
through  these  points.  The  nodes  at  the  ends  of  the  segment  are  numbered,  and 
the  segment  number  is  appended  to  the  map  near  the  midpoint  of  the  segment. 

Four  types  of  computation  are  used  by  subroutines  SHAPCOM  and  COORD  to 
produce  the  coordinates  of  points  on  a segment.  The  simplest  computation  is 
performed  for  straight  segments  and  involves  a linear  interpolation  between 
the  initial  and  final  nodes.  Another  calculation  processes  both  circular-arc 
and  S-curve  segments;  an  S-curve  is  treated  as  two  consecutive  circular  arcs. 
The  center  and  radius  of  the  circular  arc  are  found  in  subroutine  SHAPCOM. 
Subroutine  COORD  uses  the  distance  along  the  perimeter  to  compute  the  angle 
subtended  by  the  arc;  it  then  computes  the  coordinates  of  the  point  on  a ref- 
erence circle,  translates  the  reference  point,  and  rotates  the  point  to  ob- 
tain its  coordinates  on  the  segment. 

In  SHAPCOM  the  radius  of  the  circular  segment  is  found  using  a curve  fit 
to  the  solution  of  a transcendental  equation  containing  the  radius.  The  geom- 
etry for  the  derivation  of  the  equation  is  shown  in  Figure  2.  T represents 
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Figure  2.  Geometry  for  Circular  Arc  Radius  Calculation 

the  total  length  of  the  circular  segment  and  corresponds  to  FORTRAN  variable 
TOTLEN.  D is  the  length  of  a chord  connecting  the  end  points  of  the  circular 
arc  segment  and  corresponds  to  FORTRAN  variable  D in  SHAPCOM.  <p  is  one-half 
the  angle  subtended  by  the  circular  arc  at  the  center  of  the  circle.  R is 
the  radius  of  the  circular  segment.  From  fundamental  definitions  of  $ and 
si ncp , one  obtains  the  equations 

<p  - and  sin<J>  - 2p 
Dividing  these  two  equations  yields 

D si  n<t> 

T " <J> 

Since  this  equation  decreases  monotonical ly  as  <J>  increases  from  0 to  pi  (its 
maximum  permissable  value),  <p  can  be  found  by  curve  fitting  <j>  as  a function 
of  D/T.  The  curve  fit  actually  used  fits  4>2  against  1 - D/T  and  is  motivated 
by  the  approximate  solution  for  small  values  of  <j>:  replacing  sin4>/<}>  by  the 
approximation  1 - 2 / 3 ! yields  D/T  = 1 - 2 / 3 ! ; therefore  4>?  = 6(1  - D/T). 
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Values  were  tabulated  for  (<p/v)''  and  1 - sin (which  is  1 - D/T)  for 
101  values  of  <p/i r from  0 to  1 . The  relation  was  treated  as  an  odd  function, 
giving  100  more  points  for  negative  values  of  ($/tt)2,  and  a seventh-degree 
polynomial  curve  fit  was  performed.  Treating  the  relation  as  an  odd  function 
caused  all  coefficients  of  even  powers  of  (1  - D/T)  to  vanish.  Since  <p  = 

T/2R,  the  fit  for  <£  yielded  1/R.  The  result  of  the  curve  fit  is  used  by 
SHAPCOM  to  evaluate  1/R,  which  is  FORTRAN  variable  RPR. 

The  error  in  1/R  is  reduced  by  reducing  the  error  in  the  difference  of 
sin4>  = sinT/2R  and  sin<J)  = D/2R.  For  the  present  value  of  1/R  an  error  = 
sinT/2R  - D/2R  is  evaluated.  If  |e  | >0.00001,  the  following  iterative  im- 
provement is  performed.  The  error  is  positive  when  1/R  is  too  small.  A 
change  Al/R  is  selected  such  that  Al/R  = 0.0002  * sign  (e^,  and  another  error 
c2  = sin[T/2(l/R  + Al/R)]  - D/2(l/R  + Al/R)  is  evaluated.  This  error  is  either 
closer  to  zero  or  opposite  in  sign  from  e . A linear  interpolation  is  per- 
formed to  improve  the  value  of  1/R:  1/Rnew  = l/Rold  ~ £iAl/Rold^ei  ' e2)- 
When  |ej  <0.00001,  the  process  is  stopped. 

Once  the  radius  is  known,  the  center  of  the  circle  must  be  found.  There 
are  four  cases:  right  or  left  circular  arcs,  each  either  more  or  less  than 
half  a circle.  The  geometry  used  to  find  the  center  for  both  cases  of  right 
circular  arcs  is  shown  in  Figure  3.  In  this  figure,  0 is  the  slope  of  the 
line  from  starting  to  ending  point,  and  a is  +1  for  a right  circular  arc  or  -1 
for  a left  circular  arc.  H is  found  using  the  Pythagorean  theorem  on  a tri- 
angle having  sides  R and  H. 

The  computation  for  the  coordinates  (Xs,  Yg)  on  a circular  arc  at  a perim- 
eter S from  the  starting  point  (X,,  Yj)  can  be  derived  easily  by  using  vectors. 


(XCTR.YCTR) 


XCTR  = 0.5  (XNI+XE)  - aH  sine 
YCTR  = 0.5  (YNI+YE)  + aH  cose 
Note:  a = -1  for  left  circular  arcs. 


Figure  3.  Geometry  for  Circular  Arc  Center  Calculation 


is  the  vector  from  the  origin  to  the  center  of  the  circle.  The  radius  vector 


is  from  the  center  of  the  circle  to  the  starting  point.  The  rotation  matrix 

C°S  | -sin  I 

s s 

sin  p cos  p 

the  starting  point  to  the  point  on  the  arc  a 

A". 

W 

to  this  point,  and  therefore  the  coordinates  (X$,  Y ),  are  given  by 


rotates  the  radius  vector  from 
perimeter  S away.  The  vector 
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Figure  4 illustrates  the  geometry. 


To  calculate  the  coordinates  of  a point  on  a rectangular  segment,  the 
slope  components  of  the  first  side  are  determined,  and  then  appropriate  multi- 
ples of  these  components  are  added  to  the  starting  or  ending  node's  coordin- 
ates. The  distances  to  the  break  points,  BR1  and  BR2,  and  the  slope  components 
SX  and  SY  are  computed  in  subroutine  SHAPCOM.  The  geometry  for  the  calculation 
is  shown  in  Figure  5.  Using  the  FORTRAN  variable  names,  the  calculations  and 
units  are 

BR1  =0.5  (TOTLEN  - D)  miles 
BR2  = 0.5  (TOTLEN  + D)  miles 
SX  = sine  • SCR/AVMD  MCU/mile 
SY  = cose  • SCR/AVMD  MCU/mile 


where  TOTLEN  is  the  total  length  of  the  segment,  SCR  is  the  ratio  of  current 
to  overall  map  distance  conversion,  and  AVMD  is  the  current  map  distance  con- 
version. The  slope  components  and  distances  to  breaks  are  used  in  subroutine 
COORD  to  obtain  the  coordinates  (XX,  YY)  of  the  point  a perimeter  S from  the 
starting  node  of  the  rectangle.  The  following  equations  are  used: 

S < 0.95-BR1 
0.95-BR1  < S < 1.05-BR1 

1.05- BR1  < S < 0.95-BR2 
0.95-BR2  < S < 1.05-BR2 

1.05- BR2  < S 

S < 0.95-BR1 
0.95-BR1  < S < 1.05-BR1 

1.05- BR1  < S < 0.95-BR2 
0.95-BR2  < S < 1.05-BR2 

1.05- BR2  < S 


YNI-o*SY*S, 

YNI-a*  SY*  BR1 , 

YY  = YNI-a*SY*BRl+SX* (S-BR1 ) , 

YNI-a* SY* BR1+SX* (BR2-BR1 ) , 
YNF-a*SY* (BR1+BR2-S) , 


XNI+a*SX*S, 

XNI+a-SX-BRl, 

XX  = XNI+a*SX*BRl+SY* (S-BR1 ) , 
XNI+a*SX*BRl+SY* (BR2-BR1 ) , 
XNF+a*SX* (BR1+BR2-S) , 
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Figure  5.  Geometry  for  Calculating  Rectangular  Segment  Parameters 


1 


Figure  6.  Geometry  for  Calculating  Angle  Segment  Parameters 


where  (XNI,  YNI)  are  the  coordinates  of  the  starting  node,  (XNF,  YNF)  are  the 
coordinates  of  the  ending  node,  a = +1  for  a right  rectangle,  and  a = -1  for  a 
left  rectangle. 

The  coordinates  of  a point  on  an  angle  segment  are  found  by  linear  inter- 
polation between  one  end  and  the  vertex.  The  coordinates  of  the  vertex  are 

found  in  subroutine  SHAPCOM  and  are  stored  in  variables  XCTR  and  YCTR.  The 

geometry  is  shown  in  Figure  6.  Here  BR1  and  BR2  are  the  lengths  of  the  sides 

of  the  angle.  The  fraction  F of  D that  forms  a right  triangle  with  side  H and 
hypotenuse  BR1  is  found  by  eliminating  H2  from  the  two  Pythogorean  relations 
for  the  two  right  triangles  having  H as  a common  side,  and  then  solving  for  F. 
H can  then  be  found  by  using  the  right  triangle  with  hypotenuse  BR1 . The  re- 
sults are 

F = [1  - ( BR22 -BR1 2 ) /D 2 ] / 2 


H = a VbRI  2-  (F*D) 2 


where  n = +1  for  an  angle  that  lies  to  the  right  of  the  line  connecting  the 
end  points,  or  a = -1  for  an  angle  on  the  left  side  of  that  line. 

The  coordinates  of  the  vertex  of  the  angle  are 

XCTR  = XNI+(cosG-F-D-sinQ»H)*SCR/AVMD 
YCTR  = YN 1+ ( s i nG • F • D- cos0 • H ) • SCR/AVMD 

where  SCR,  AVMD,  XNI,  and  YNI  have  the  same  meanings  as  for  a rectangular  seg- 
ment. The  coordinates  (XX,  YY)  of  a point  S miles  from  the  starting  node  of 
the  angles  are  found  in  subroutine  COORD  by  linear  interpolation.  The  equa- 
tions used  are  equivalent  to 


XN I + ( XCTR-XN I ) • S/ BR 1 , 

s < 

BR1 

_XCTR+ (XNF-XCTR ) • (S-BR1  )/BR2, 

BR1 

< S 

~YNI+(YCTR-YNI)-S/BR1, 

s < 

BR1 

_YCTR+(YNF-YCTR) • ( S- BR  1 )/BR2 , 

BR1 

< S 

2.  DATA  STORAGE 

Three  files  generated  by  program  RCINPT  are  saved  on  disk  for  use  by  later 
programs:  segment  data  on  file  TAPE1 , node  data  on  file  TAPE2,  and  street  data 
on  file  TAPE3.  Files  TAPE1  and  TAPE2  are  used  by  programs  PHASE2,  PHASE3,  and 
PHASE4.  File  TAPE3  is  used  only  by  program  PHASE4. 

The  segment  data  are  stored  in  array  STG,  which  is  in  blank  COMMON.  The 
array  is  equivalenced  with  array  ISTG  so  that  both  integer  and  floating-point 
segment  data  can  be  accessed  by  reference  to  the  appropriate  array  name.  All 
of  the  data  in  arrays  STG  and  ISTG  come  from  the  input  map-description  records. 
In  the  following,  the  variables  NSTR , NN1 , NN2,  LEN,  NH,  NSPD,  NWAY,  NRQF, 
NXMID,  NYMID,  and  NSF  represent  numeric  subscripts  1 through  11,  respectively. 

The  street  number  and  starting  node  of  each  string  are  read  by  the  call 
to  LINEIN  at  statement  number  210  in  main  program  RCINPT.  The  initial  node 
number  is  stored  in  ISTG(NNl.KI)  two  statements  after  statement  number  240. 
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The  street  number  is  saved  in  variable  NUMST  until  later.  KI  is  the  line  num- 
ber at  the  start  of  string  processing.  Each  segment  is  read  by  LINEIN  at 
statement  250.  The  length  and  the  number  of  houses  on  the  right  side  are 
stored  in  STG(LEN,KF)  and  ISTG(NH.KF)  after  statement  260.  KF  is  the  current 
line  number.  The  number  of  houses  on  the  left  side  is  saved  in  array  NHL  for 
later  use.  The  next  node  number  is  stored  in  both  ISTG(NN2,KF)  and 
ISTG(NN1 ,KF+1 ) , six  statements  after  number  260.  If  the  string  continues,  the 
node  is  both  the  ending  node  of  one  segment  and  the  starting  node  of  the  next. 
If  the  string  ends  on  the  node,  the  next  string  will  overwrite  ISTG(NN1 ,KF+1 ) 
with  its  starting  node  number.  The  speed  limit,  number  of  ways  of  travel,  num- 
ber of  sides  serviced  on  one  pass,  and  refuse  quantity  adjustment  factor  are 
read  by  the  call  to  LINEIN  two  statements  before  statement  280.  These  param- 
eters, or  default  values  if  the  parameters  are  missing,  apply  to  all  segments 
in  the  string.  The  speed  limit,  number  of  ways  of  travel,  and  refuse  quantity 
adjustment  factor  are  stored  in  STG(NSPD,K),  ISTG(NWAY ,K) , and  STG(NRQF,K)  in 
the  loop  through  statement  290.  K is  the  loop  index  and  has  values  equal  to 
the  line  numbers  of  the  segments.  Also  in  this  loop,  the  street  number  is 
stored  in  ISTG(NSTR,K) ; and  if  both  sides  of  the  street  are  serviced  on  one 

pass,  the  houses  on  the  left  side  are  added  to  those  on  the  right  side. 

/ 

The  midpoint  coordinates  of  the  segments  and  the  segment  shape  codes  are 
stored  in  STG(NXMID,K) , STG(NYMID,K) , and  ISTG(NSF.K)  in  the  loop  through 
statement  690  as  the  segments  are  plotted.  If  collection  is  from  only  one 
side  of  a street  at  a time,  additional  segments  are  generated  in  the  loop 
through  statement  700.  These  additional  segments  are  made  one-way  segments  if 
the  corresponding  original  segment  was  a two-way  segment.  The  signs  on  the 
number  of  houses  on  the  original  and  new  segments  are  made  minus  to  indicate 
collection  from  one  side  of  the  segment  at  a time. 

When  all  plotting  is  complete,  the  count  of  segments  KF,  as  much  of  the 
ISTG  array  as  has  been  filled,  and  the  map  distance  conversion  factor  for  the 
first  input  map  SVAV(l)  are  written  to  TAPE1  by  the  binary  WRITE  following 
statement  1030.  SVAV(l)  is  set  at  statement  90  and  is  computed  from  a product 
using  the  map  scale  and  the  coordinate  scale,  both  of  which  are  read  from  the 
first  card  in  the  map-description  record. 


The  data  written  to  TAPE2  consist  of  variables  NHTOT,  TOTREF,  and  KNODES, 
and  node  data  arrays  NODNUM,  NBS,  XNOD,  and  YNOD.  Variable  KNODES  and  the  four 
arrays  are  stored  in  COMMON  block  NDDATA.  The  total  number  of  houses  (NHTOT) 
and  the  total  refuse  (TOTREF)  are  accumulated  in  the  loop  through  statement  290, 
using  data  from  the  map-description  strings.  KNODES,  the  count  of  nodes,  is  in- 
cremented at  the  statements  before  statements  245  and  265,  immediately  follow- 
ing the  insertion  of  node  numbers  into  array  NODNUM.  This  array  is  kept  in  in- 
creasing order  by  using  function  IFIND  to  determine  where  the  node  number 
should  go  and  subroutine  M0VE5  to  open  a space  for  it  in  the  NODNUM  array. 

NBS(I)  contains  up  to  six  numbers  of  segments  bounding  the  node  NODNUM(I).  The 
segment  number  is  appended  to  that  entry  in  array  NBS  corresponding  to  the 
starting  node  following  statement  260  and  corresponding  to  the  ending  node  at 
statement  265.  The  x-  and  y-coordi nates  of  the  node  are  stored  in  XNOD  and 
YNOD  near  the  end  of  the  loop  through  statement  690.  The  binary  WRITE,  two 
statements  after  1030,  writes  NHTOT,  TOTREF,  KNODES,  and  the  filled  parts  of 
arrays  NODNUM,  NBS,  XNOD,  and  YNOD  to  file  TAPE2. 

The  street  numbers  and  names  are  read  from  the  first  record  of  card  in- 
put by  subroutine  STRINP.  These  data  are  stored  in  arrays  NUMSTR  and  NAMSTR 
in  blank  COMMON.  After  100  street  name  cards  are  read,  the  NUMSTR  and  NAMSTR 
arrays  are  written  to  TAPE3  by  a BUFFER  OUT  statement  in  subroutine  STRINP. 

When  the  last  street  name  card  is  read,  all  of  the  NUMSTR  and  NAMSTR  arrays 
are  buffered  out  to  TAPE3,  so  the  records  on  TAPE3  are  always  the  same  length. 
The  unused  portions  of  these  arrays  will  be  filled  with  zeros.  These  arrays 
are  not  needed  after  STRINP  returns  control  to  the  main  program.  Since  NUMSTR 
and  NAMSTR  occupy  the  same  storage  as  the  STG  array,  they  will  be  overwritten 
with  segment  data  after  control  returns  to  the  main  program. 

3.  PURPOSE  AND  PERFORMANCE 


In  this  section  the  simplest  subroutines  are  described  first  so  their 
workings  will  be  clear  when  they  are  mentioned  again  in  the  descriptions  of 
the  more  complicated  subroutines  and,  finally,  of  the  main  program.  Logic 
flowcharts  are  given  in  Appendix  A.  Complete  program  listings  are  provided 
in  Appendix  B.  In  Appendix  C,  the  more  important  variables  in  the  following 
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subroutines  are  defined  in  terms  of  their  specific  meaning  for  each 
subroutine. 

a.  Subroutine  M0VE5 

The  purpose  of  subroutine  M0VE5  is  to  move  all  data  starting  at  and 
following  a given  subscript  to  a position  in  the  array  starting  at  a given 
higher  subscript.  Argument  II  represents  the  subscript  from  which  the  data 
will  be  moved,  and  IF  represents  the  subscript  to  which  the  data  will  be  moved. 
Arguments  A1 , A2,  A3,  A4,  and  A5  are  the  names  of  the  arrays  in  which  data  will 
be  moved.  Subroutine  M0VE5  starts  by  testing  whether  the  data  are  to  be  moved 
to  a higher  subscript.  If  the  final  subscript,  IF,  is  less  than  the  initial 
subscript,  II,  the  program  returns  control  to  the  calling  program.  If  the 
initial  and  final  subscripts  are  in  proper  relation,  the  subroutine  executes 
the  loop  through  statement  10,  which  moves  each  item,  starting  at  subscript  II, 
the  distance  necessary  to  reach  IF.  All  data  between  subscripts  II  and  IF  are 
moved,  beginning  with  the  last  item.  Motion  starts  at  the  higher  subscript  so 
that  no  data  that  must  be  moved  later  are  overwritten.  Zeros  are  stored  at 
the  location  corresponding  to  subscript  II,  from  which  data  were  moved.  Con- 
trol is  then  returned  to  the  calling  program. 

b.  Function  IFIND 

Function  IFIND  uses  a binary  search  to  locate  a given  number  in  an 
array;  the  subscript  corresponding  to  the  location  of  the  number  is  assigned 
as  the  value  of  IFIND.  If  the  number  is  not  found,  the  function  sets  a value 
for  IFIND  equal  to  the  negative  of  the  subscript  at  which  the  number,  to  be  in 
numerical  order,  should  be  inserted.  (The  array  is  assumed  to  be  in  increasing 
order.)  The  comment  cards  at  the  beginning  of  function  IFIND  list  the  latest 
changes  to  the  function  and  state  the  function's  purpose. 

Argument  NUM  is  the  number  that  is  sought  in  array  I ARRAY.  The 
length  of  array  IARRAY  is  given  by  argument  LEN.  Function  IFIND  begins  by 
checking  that  LEN  > 0.  If  LEN  < 0,  the  function  assigns  a value  for  IFIND  of 
-1.  This  value  indicates  that  the  number  sought  is  not  in  the  array  and  would 
be  stored  as  the  first  entry  in  the  array.  The  binary  search  uses  variables 
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II,  IP,  and  IF  as  pointers.  II  is  the  subscript  of  the  front  of  the  region 
being  searched,  IP  is  the  subscript  of  the  item  being  compared  to  the  number 
sought,  and  IF  is  the  subscript  of  the  last  item  in  the  region  being  searched. 
Variable  II  is  initially  set  to  1 at  statement  5,  and  variable  IF  is  set  to 
the  end  of  the  array  in  the  next  statement.  The  pointer  IP  is  the  subscript 
about  midway  between  II  and  IF. 

The  computation  of  IP  occurs  at  statement  10.  The  statement  follow- 
ing statement  10  compares  the  number  being  sought,  NUM,  to  the  data  at 
I ARRAY (IP).  If  NUM  < IARRAY(IP),  control  transfers  to  statement  20,  indica- 
ting that  the  number  is  in  the  front  half  of  the  region  being  searched;  at 
statement  20  the  final  pointer  is  moved  to  the  subscript  preceeding  the  point 
just  searched.  If  NUM  > IARRAY(IP),  control  transfers  to  statement  30,  indi- 
cating that  the  number  being  sought  follows  the  subscript  just  inspected;  at 
statement  30  the  initial  pointer,  II,  is  set  to  the  present  pointer,  IP,  plus 
1.  If  the  number  sought  is  found  at  IARRAY(IP),  control  transfers  to  state- 
ment 50,  where  IFIND  is  set  equal  to  the  current  pointer  and  control  returns 
to  the  calling  program.  Where  NUM  is  unequal  to  I ARRAY (IP) , control  resumes 
at  statement  40  after  the  initial  or  final  pointers  are  moved.  At  statement 
40  the  final  pointer  is  compared  to  the  initial  pointer;  if  IF  > II,  control 
is  transferred  to  statement  10. 

At  statement  10  the  search  is  resumed  on  the  appropriate  half  of  the 
region  examined  previously.  If  the  final  pointer  becomes  less  than  the  initial 
pointer,  the  number  sought  is  not  in  the  table.  In  this  case,  control  resumes 
following  statement  40,  and  the  value  of  IFIND  is  set  to  the  negative  of  the 
current  pointer.  If  the  number  at  the  current  pointer  is  less  than  the  number 
being  sought,  IFIND  is  set  to  -(IP  + 1 ) so  the  number  can  be  inserted  in  the 
appropriate  place.  Control  then  returns  to  the  calling  program. 

c.  Subroutine  STRINP 

Subroutine  STRINP  reads  and  prints  the  first  record  of  data:  the 
title  card  and  the  street-name  information.  It  also  writes  the  street-name 
data  on  file  TAPE3.  It  has  one  argument,  NIIR,  which  indicates  to  the  call- 
ing program  the  presence  or  absence  of  street-name  data. 
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Blank  COMMON  is  used  to  hold  the  title  and  the  street  numbers  and 
names.  The  title  is  left  permanently  in  blank  COMMON,  but  the  street  numbers 
and  names  are  overwritten  later.  Variable  MSTINC  is  set  to  100  in  a DATA 
statement.  This  variable  controls  the  maximum  number  of  streets  stored  in 
core  before  the  data  are  written  to  TAPE3. 

The  subroutine  initially  assumes  that  street-name  data  will  be  pres- 
ent and  sets  NIIR  = 2 to  indicate  this.  Variable  NS  is  set  to  0 and  is  used 
to  count  the  number  of  cards  read.  The  title  card  is  read  according  to  the 
format  in  statement  5:  8A10.  The  storage  set  aside  for  the  street  numbers  is 

set  to  0 by  the  loop  through  statement  20.  The  next  statement  attempts  to 
read  100  cards,  each  containing  a street  number  and  a street  name.  If  an  end- 
of-record  card  is  encountered  before  100  cards  have  been  read,  the  zeros  stored 
in  the  street-number  array  will  remain  unchanged  after  the  end-of-record  is  en- 
countered. 

The  loop  through  statement  40  starts  at  the  end  of  the  storage  set 
aside  for  street  numbers  and  searches  for  a nonzero  value.  When  a nonzero 
street  number  is  encountered,  control  transfers  to  statement  50.  (Variable 
INC  is  a count  of  the  number  of  cards  read.)  If  this  loop  finds  no  nonzero 
street  numbers,  control  transfers  to  statement  120. 

At  statement  50  a loop  through  statement  80  prints  the  number  and 
name  of  each  street  read  in.  After  statement  80,  the  information  in  arrays 
NIJMSTR  and  NAMSTR  (street  numbers  and  names)  is  written  to  TAPE3  by  a BUFFER 
OUT  statement  using  parity  1.  The  unit  status  is  tested  by  the  next  state- 
ment; when  a parity  error  is  encountered,  control  is  transferred  to  statement 
90,  which  prints  an  error  message.  In  either  case  control  then  resumes  in 
statement  110,  which  increments  the  cumulative  count  of  streets,  NS,  by  the 
number  of  streets  read  and  written. 

Statement  120  tests  for  an  end-of-record  in  the  card  input.  If  no 
end-of-record  is  encountered,  control  transfers  to  statement  10,  and  additional 
cards  are  read.  If  an  end-of-record  is  encountered,  control  resumes  at  state- 
ment 130.  If  no  streets  have  been  read  by  this  point,  an  error  message  is 
printed  and  variable  NIIR  is  set  to  1 . If  streets  are  found,  control  transfers 
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from  statement  130  to  statement  150.  After  statement  150  an  end-of-file  is 
written  on  TAPE3  and  the  file  is  rewound.  Control  then  returns  to  the  calling 
program. 


d.  Subroutine  NUMBER 

Subroutine  NUMBER  appends  numbers  to  plotted  output.  Its  use  is  al- 
most identical  to  that  of  the  standard  Calcomp  number  routine,  the  primary  dif- 
ference being  that  the  last  argument  in  subroutine  NUMBER  gives  an  alphanumeric 
format  rather  than  an  integer  format  code. 

Subroutine  NUMBER  has  six  arguments.  The  first  two  give  the  coor- 
dinates, in  plotter  inches,  of  the  lower  left  corner  of  the  field.  The  third 
gives  the  height,  in  inches,  of  the  digits.  The  fourth  is  the  number  to  be 
plotted.  The  fifth  is  the  angle  at  which  the  number  is  to  be  plotted,  meas- 
ured in  degrees  counterclockwise  from  the  horizontal.  The  last  argument  is  an 
alphanumeric  format  up  to  10  characters  long,  which  describes  the  appearance 
of  the  plotted  number. 

Array  TEXT  is  used  to  hold  the  character  representation  of  the  num- 
ber. Up  to  30  characters  are  allowed.  The  first  executable  FORTRAN  statement 
sets  this  array  to  three  words  of  blanks.  The  second  statement  moves  the  for- 
mat into  the  second  word  of  the  array  FORM.  The  first  and  third  words  of  this 
array  have  been  preset  to  a left  and  a right  parenthesis  by  a DATA  statement. 
The  ENCODE  statement  converts  the  number  from  binary  form  in  variable  NUM  to 
character  form  in  array  TEXT,  according  to  format  FORM.  A character  count, 
variable  NC,  is  set  to  30. 

The  loop  through  statement  10  searches  for  the  last  nonblank  char- 
acter in  array  TEXT.  Each  time  a blank  is  found,  starting  at  the  end  of  the 
TEXT  array,  the  character  count  (NC)  is  decremented  by  1.  When  a nonblank 
character  is  encountered,  control  transfers  to  statement  20 

Statement  20  calls  the  standard  SYMBOL  subroutine  to  plot  the  char- 
acter representation  of  the  number.  Control  then  returns  to  the  calling  pro- 
gram. 


I: 
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e.  Subroutine  AXIS 

Subroutine  AXIS  draws  an  axis  with  tic  marks,  numbers,  and  a label. 

The  axis  and  numbers  may  be  drawn  at  any  specified  angle,  and  the  format  for 
the  numbers  must  be  specified.  The  axis  can  be  drawn  with  tic  marks  extending 
through  it  or  with  ha  If- size  tic  marks  on  the  numbered  side  of  the  axis.  An- 
other option  allows  the  tic  marks,  alone,  to  be  plotted  (without  the  axis). 

Subroutine  AXIS  has  13  arguments.  The  first  two  arguments,  X and  Y, 
are  the  x-  and  y-coordinates  of  the  low-value  end  of  the  axis,  in  plotter 
inches.  The  next  two  arguments,  VI  and  VF,  are  the  values  of  the  start  and 
end  of  the  axis.  The  fifth  argument,  SCALE,  is  the  scale  of  the  axis  in  plot- 
ter inches  per  axis  unit.  The  sixth  argument,  TIC,  is  the  value  of  the  inter- 
val between  small  tic  marks.  The  seventh  argument,  DLBL,  is  the  value  of  the 
interval  between  the  larger  numbered  tic  marks. 

The  eighth  argument,  MODE,  is  the  mode  of  plotting  of  the  axis.  If 
MODE  = 0,  the  axis  is  plotted  with  tic  marks  extending  across  the  axis.  If 
MODE  = 1,  only  tic  marks  are  plotted.  If  MODE  = 2,  the  axis  is  plotted  with 
tic  marks  on  only  the  numbered  side  of  the  axis.  If  MODE  = 3,  half-size  tic 
marks  are  plotted.  Multiples  of  4 can  be  added  to  the  mode  to  decrease  the 
number  of  numbered  tic  marks.  If  0 or  4 is  added  to  the  mode,  every  large 
tic  mark  corresponding  to  an  interval  of  DLBL  will  be  numbered.  If  8 is  added 
to  the  mode,  every  second  large  tic  mark  will  be  numbered.  If  12  is  added  to 
the  mode,  every  third  large  tic  mark  is  numbered. 

The  ninth  argument,  FMT,  is  the  format  of  the  numbering  on  the  axis. 
The  tenth  argument,  ANGAX,  is  the  angle  of  the  axis  measured  in  degrees  counter- 
clockwise from  the  right  horizontal  axis.  For  ANGAX  < 0,  the  label  and  numbers 
go  on  the  clockwise  side  of  the  axis.  For  ANGAX  > 0,  the  label  and  numbers  go 
on  the  counterclockwise  side  of  the  axis.  The  label  runs  in  the  same  direction 
as  the  axis  if  j ANGAX | < 360°,  but  in  the  opposite  direction  for  360°  < 

| ANGAX | < 720°. 

The  eleventh  argument,  ANGNM,  is  the  angle  of  the  numbers,  in  de- 
grees. The  twelfth  argument,  LBL,  is  an  array  containing  the  label  for  the 
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axis.  The  thirteenth  argument,  NC,  is  the  number  of  characters  in  array  LBL 
to  be  used  in  the  label.  NC  may  equal  0,  in  which  case  no  label  is  appended 
to  the  axis. 

The  first  executable  statement  in  AXIS  sets  variable  IFP  to  1,  assum- 
ing that  the  format  specifies  floating-point  numbers.  If  this  assumption  is 
incorrect,  IFP  will  be  set  to  0 when  the  format  is  scanned.  The  cosine  and 
sine  of  the  axis  angle  are  computed  next.  If  either  is  within  0.0001  of  0, 
it  is  set  to  0.  If  either  is  close  to  1 or  -1 , it  is  set  to  the  appropriate 
value.  Thus  roundoff  errors  are  eliminated  from  axes  that  occur  at  some  mul- 
tiple of  90°. 

The  next  statement  evaluates  TCC  and  TCS,  which  are  the  horizontal 
and  vertical  components  of  the  small  tic  marks.  The  next  statement  evaluates 
variable  IPEN.  This  variable  will  be  equal  to  2 if  the  entire  axis  is  to  be 
drawn,  or  equal  to  3 if  only  the  tic  marks  are  to  be  drawn.  Logical  variables 
CCS  and  CS  are  set  to  TRUE,  assuming  that  the  tic  marks  will  be  drawn  on  both 
the  counterclockwise  and  clockwise  sides  of  the  axis.  The  next  statement  tests 
the  two  bit  of  the  mode  to  see  whether  the  tic  marks  extend  through  the  axis. 

If  the  two  bit  of  the  mode  is  0,  control  transfers  to  statement  10.  If  the 
two  bit  is  1,  variable  CS  is  set  to  TRUE  or  FALSE  according  to  whether  or  not 
the  numbering  is  on  the  clockwise  side  of  the  axis.  CCS  is  set  to  the  comple- 
ment of  CS. 

Statement  10  sets  variable  N to  the  total  number  of  tic  marks  to  be 
drawn  on  the  axis.  The  next  statement  positions  the  plotter  pen  at  the  begin- 
ning of  the  axis.  Variable  M is  set  to  the  number  of  tic  marks  between  one 
large  tic  mark  and  the  next.  The  next  FORTRAN  line  sets  variables  XX  and  VY 
equal  to  the  initial  position  of  the  pen.  These  variables  will  contain  the 
current  pen  position  as  the  axis  is  drawn.  Variables  XINC  and  YINC  are  com- 
puted; these  are  the  horizontal  and  vertical  increments  during  the  drawing  of 
the  axis. 


The  loop  through  statement  20  draws  the  axis.  Each  pass  through  the 
loop  draws  the  axis  from  one  tic  mark  to  the  next.  Variable  F indicates 
whether  the  tic  mark  is  a small  or  a large  one.  For  each  small  tic  mark,  F = 1; 
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F = 2 for  each  large  tic  mark.  The  pen  is  positioned  at  coordinates  XX  and  YY 
If  variable  CCS  is  TRUE,  a tic  mark  is  drawn  on  the  counterclockwise  side  of 
the  axis.  If  variable  CS  is  TRUE,  a tic  mark  is  drawn  on  the  clockwise  side 
of  the  axis.  The  next  statement  returns  the  pen  to  coordinates  XX  and  YY, 
which  are  on  the  axis.  The  last  two  statements  in  the  loop  increment  coordin- 
ates XX  and  YY . The  call  to  PLOT  following  statement  20  advances  the  pen  to 
the  high-value  end  of  the  axis.  The  next  call  to  PLOT  returns  the  pen  to  the 
beginning  of  the  axis,  causing  the  central  line  of  the  axis  to  be  drawn  twice 
unless  only  tic  marks  are  to  be  drawn.  The  next  statement  tests  variable  IPEN 
to  see  whether  only  tic  marks  are  to  be  drawn.  In  this  case,  control  returns 
to  the  calling  program.  If  the  entire  axis  is  drawn,  numbers  will  now  be  ap- 
pended to  the  large  tic  marks. 

Variable  DIFA  is  the  difference,  measured  in  radians,  between  the 
angle  of  the  numbers  and  the  angle  of  the  axis.  The  next  statement  computes 
the  sine  and  cosine  of  this  angle.  The  cosine  of  twice  the  angle  is  also  com- 
puted. The  next  line  sets  variable  S to  either  1 or  -1 , depending  on  whether 
or  not  the  axis  angle  is  positive.  Variable  NCN,  the  number  of  characters  in 
the  numbering,  is  initially  set  to  0. 

The  loop  through  statement  22  will  scan  the  format  for  the  output 
type.  If  the  output  is  not  of  type  E,  F,  G,  or  I,  the  program  stops  with  num- 
bered STOP  567.  If  an  I-type  number  is  specified,  control  transfers  to  state- 
ment 23,  where  variable  IFP  is  set  to  0.  If  type  £,  F,  or  G is  found,  control 
transfers  to  statement  24. 


At  statement  24  a loop  is  started  through  statement  26,  which  scans 
the  remaining  characters  to  find  the  field  width.  The  field  width  is  stored 
in  variable  NCN.  When  the  scan  finds  a noninteger  character,  control  trans- 
fers to  statement  28. 

Following  statement  28,  variable  HWD.  the  half-width  of  the  number, 
is  evaluated.  The  position  of  the  beginning  of  the  number  field  relative  to 


the  corresponding  point  on  the  axis  is  determined  next  in  terms  of  the  compon- 
ents normal  and  tangential  to  the  axis  (DIFN  and  DIFT).  The  x-  and  y-components 
of  this  position  difference  are  evaluated  next  in  variables  DELX  and  DELY.  The 


spacing  between  large  tic  marks  is  computed  in  variables  XINC  and  YINC.  Vari- 
able N is  the  number  of  large  tic  marks.  Variable  I INC  is  computed  next  and 
gives  the  number  of  large  tic  marks  between  numbered  large  tic  marks. 

The  loop  through  statement  30  appends  numbers  to  the  appropriate 
large  tic  marks.  The  value  of  the  number  is  kept  in  variable  V;  if  the  num- 
ber is  to  be  an  integer,  the  value  is  changed  to  integer  preceeding  state- 
ment 30.  The  number  is  appended  to  the  plot  at  statement  30.  Variables  V and 
IV  are  equivalenced  so  that  the  integer  value  can  be  accessed  as  variable  V. 

The  statement  following  statement  30  tests  to  see  whether  a blank 
label  or  no  label  is  to  be  appended  to  the  axis.  In  either  case,  control  re- 
turns to  the  calling  program.  Otherwise,  the  remaining  statements  determine 
the  position  and  direction  of  the  label  for  the  axis.  Variable  S is  set  to 
1,  indicating  that  the  label  runs  in  the  same  direction  as  the  axis.  The 
next  statement  changes  S to  -1  if  the  angle  of  the  axis  is  an  odd  number  of 
full  rotations  from  the  interval  -360.0001  to  360.0001.  The  coordinates  XX 
and  YY  of  the  midpoint  of  the  label  are  then  computed.  If  the  axis  angle  is 
0,  the  y-coordinate  of  the  midpoint  is  recomputed.  The  beginning  coordinate 
of  the  lower  left  corner  of  the  label  field  is  computed  next  in  variables  XXX 
and  YYY.  The  label  is  then  plotted  in  characters  0.15  inch  high  by  the  call 
SYMBOL.  Control  then  returns  to  the  calling  program. 

f.  Subroutine  MAPGRID 

Subroutine  MAPGRID  plots  a grid  around  the  map  of  the  collection  re- 
gion so  that  coordinates  may  be  easily  read  from  the  map.  Arguments  XMIN  and 
XMAX  are  the  minimum  and  maximum  coordinates  of  the  map  region  to  be  plotted. 
Argument  XLEN  is  the  length  of  the  x-axis,  in  inches.  Arguments  YMIN,  YMAX, 
and  YLEN  are  the  corresponding  y-direction  parameters.  If  a drum  plotter  is 
used,  the  last  argument,  YHCUT,  is  the  height  of  strips  of  the  map.  If  a 
flatbed  plotter  is  used,  YHCUT  should  be  the  maximum  height  allowed  for  plots. 

The  first  FORTRAN  line  evaluates  XDEL  and  YDEL,  which  are  the  number 
of  map  coordinate  units  (MCU)  spanning  the  map  in  the  x-  and  y-directions . The 
next  line  evaluates  XSC  and  YSC,  which  are  the  x-  and  y-scales,  in  inches  per 
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I1CU.  Y INC , the  number  of  units  in  the  y-direction  in  one  strip  of  the  map, 
is  evaluated  next  as  the  height  of  a strip  divided  by  the  y-scale.  If  the 
height  of  a strip  is  0,  YINC  is  set  to  the  number  of  units  in  the  y- 
direction.  If  YHCUT  is  0,  the  plotter  paper  is  assumed  to  be  high  enough  to 
accommodate  the  entire  grid.  If  the  spread  in  the  y-direction  exceeds  one 
unit,  YINC  is  set  to  its  integer  part.  YHCUT  is  reset  so  that  one  strip  of 
the  map  includes  an  integral  number  of  coordinate  units  in  the  y-direction. 
This  number  is  used  by  the  other  map-plotting  subroutines. 

IDELX  and  IDELY  are  computed  next.  They  are  the  number  of  units 
between  tic  marks  drawn  in  the  vertical  and  horizontal  directions.  The 
marks  are  at  most  5 inches  apart  to  allow  coordinates  to  be  read  accurately. 
NPL  is  the  number  of  pieces  long  the  plot  will  be.  NH  is  one  greater  than 
the  number  of  y-direction  MCU  spanned  by  the  bottom  and  top  horizontal  axes 
of  each  piece.  NV  is  one  greater  than  the  number  of  x-direction  MCU 
spanned  by  the  leftmost  and  rightmost  vertical  axes  of  each  piece. 

The  loop  through  statement  50  draws  the  horizontal  and  vertical 
axes  on  each  piece  of  the  overall  grid.  XDISPL  is  the  displacement  in  the 
x-direction  to  the  beginning  of  the  Jth  piece.  YSTOP  and  Y START  are  the 
final  and  initial  values  of  the  y-axis.  The  axis  angle,  ANGAX,  is  set  to  0 
for  the  horizontal  axes.  The  80-character  problem  title  will  be  used  as  a 
label  on  the  first  horizontal  axis.  Variable  NC  is  set  for  80  characters. 

The  loop  through  statement  20  draws  horizontal  axes  for  one  piece 
of  the  map.  The  axis  mode,  M,  is  set  to  0,  indicating  that  the  full  axis  is 
to  be  drawn.  The  next  statement  transfers  control  to  statement  10  if  the 
first  or  last  axis  is  being  drawn.  The  next  statement  bypasses  axis  plot- 
ting if  the  y-coordinate  value  represented  by  1-1  is  not  a multiple  of  IDELX. 
The  axis  mode,  M,  is  set  to  1 for  all  axes  within  the  boundary  of  the  piece, 
causing  only  tic  marks  to  be  plotted. 

At  statement  10,  the  y-coordinate  of  the  beginning  of  the  axis  is 
computed.  The  axis  is  then  plotted.  After  the  first  call  to  AXIS,  NC  is 
set  to  0 so  that  the  label  will  not  be  appended  on  subsequent  calls. 
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At  statement  20  the  axis  angle  is  set  to  360°.  The  numbering  will 


now  appear  above  the  last  axis  (the  one  at  the  top  of  the  piece).  After  the 
loop  through  statement  20,  the  axis  angle  is  changed  to  90°. 

The  loop  through  statement  40  draws  vertical  axes.  The  mode,  M,  is 
initially  set  to  0.  If  the  first  or  last  vertical  axis  is  being  drawn,  con- 
trol  transfers  to  statement  30.  If  not,  the  mode  is  reset  to  1 . If  the  x- 
coordinate  represented  by  1-1  is  not  a multiple  of  IDELY,  axis  plotting  is  by- 
passed. 

At  statement  30,  the  x-position  of  the  axis  relative  to  the  beginning 
of  the  piece  is  computed.  In  the  next  statement  subroutine  AXIS  is  called,  and 
the  displacement  from  the  beginning  of  the  first  piece  to  the  beginning  of  the 
current  piece  is  added  to  XH  to  give  the  appropriate  starting  x-coordinate.  At 
statement  40  the  axis  angle  is  changed  to  -270°  so  that  the  last  vertical  axis 
will  be  numbered  on  the  right  side.  Statement  50  is  the  end  of  the  loop  that 
draws  each  piece  of  the  overall  map.  Control  returns  to  the  calling  program. 

g.  Subroutine  SHAPC0M 

Subroutine  SHAPC0M  sets  up  parameters  in  COMMON  block  COPARM  that  de- 
scribe the  geometrical  properties  of  a segment.  These  parameters  are  used  by 
subroutine  COORD  to  produce  the  coordinates  of  points  on  a segment. 

Subroutine  SHAPCOM  has  five  arguments.  Argument  TOTLEN  gives  the 
total  length  of  the  segment,  in  miles.  Argument  AVMD  gives  the  number  of 
miles  per  MCU  on  the  overall  map.  Argument  CNVLEN  gives  the  number  of  miles 
per  length  unit.  Argument  SCR  gives  a scale  ratio  for  AVMD  on  the  present 
map  compared  to  AVMD  on  the  first  map.  Argument  MODE  indicates  whether  data 
cards  should  be  printed  if  errors  are  encountered.  MODE  is  always  1 in  the 
current  version  of  the  program  because  the  data  cards  are  printed  when  read 
by  RCINPT.  The  values  of  the  arguments  are  sent  to  subroutine  SHAPCOM,  and 
all  output  values  from  SHAPCOM  are  placed  in  COMMON  block  COPARM. 
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In  COMMON  block  COPARM,  variable  SF  indicates  the  shape  of  the  seg- 
ment. XNI  and  XNF  are  the  x-coordinates  of  the  initial  and  final  nodes  of  the 
, segment.  YNI  and  YNF  are  the  y-coordi nates  of  these  nodes.  SX  and  SY  are  the 

slope,  in  MCU  per  mile,  in  the  x and  y directions.  RPR  is  the  reciprocal  of 
the  radius  of  curvature  for  circular  segments  and  the  circular  portions  of  S- 
curves.  Cll  and  C12  are  the  position  differences  in  MCU  of  the  starting  point 
and  center  of  a circular  arc  or  of  the  first  half  of  an  S-curve.  XCTR  and  YCTR 
are  the  center  coordinates  in  MCU  for  a circular  arc  or  half  an  S-curve.  BR1 
is  the  distance  in  miles  from  the  start  of  a segment  to  some  particular  point 
on  that  segment.  It  is  not  used  for  straight  segments.  For  circular  segments, 
BR1  is  the  total  perimeter.  For  an  S-curve,  BR1  is  the  perimeter  to  the  mid- 
point of  the  S-curve.  For  a rectangular  segment,  BR1  is  the  distance  to  the 
first  bend  in  the  rectangle.  For  an  angle,  BR1  is  the  distance  to  the  vertex. 
BR2  is  defined  only  for  rectangular  segments  and  angles.  It  is  the  distance 
in  miles  from  the  start  of  a rectangular  segment  to  the  second  bend.  For  an 
angle,  BR2  is  the  length  of  the  second  side.  SGN  is  -1  for  shapes  involving 
the  L (left)  prefix;  otherwise,  SGN  is  +1. 

Subroutine  SHAPCOM  begins  execution  by  assuming  that  the  shape  code 
indicates  a straight  line.  Break  indicators  BR1  and  BR2  are  set  to  0.  DX 
and  DY,  the  x-  and  y-components  of  the  vector  from  the  initial  to  the  final 
nodes  on  the  segment,  are  computed.  The  x-  and  y-components  of  the  slope  of 
the  vector,  measured  in  MCU  per  mile,  are  computed  and  stored  in  SX  and  SY. 

The  shape  code  is  tested;  if  the  segment  proves  to  be  a straight  line  or  is 
not  to  be  plotted,  the  subroutine  returns  control  to  the  calling  program.  For 
any  other  shape  code,  execution  continues.  The  angle  of  the  vector  from  the 
starting  to  the  stopping  node  is  computed  as  variable  THETA.  The  distance  from 
the  starting  to  the  stopping  point,  D,  is  computed  in  miles.  If  the  shape  code 
indicates  a shape  other  than  circular  or  S-curve,  control  transfers  to  state- 
ment 60.  If  the  distance  from  the  starting  to  the  stopping  node  is  less  than 
the  total  perimeter  of  the  segment,  control  transfers  to  statement  45.  Other- 
wise,  the  shape  code  is  set  to  0 and  control  returns  to  the  calling  program. 

At  statement  45  the  coordinates  of  the  final  node  are  stored  in  vari- 
ables XE  and  YE.  The  first  break,  BR1 , is  set  to  the  total  length  of  the  seg- 
ment. Variable  DD  is  set  to  the  straight-1 ine  distance  from  the  starting  to 
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the  stopping  node.  If  the  shape  code  indicates  a circular  segment,  control 
transfers  to  statement  50.  If  not,  variables  XE  and  YE  are  reset  to  the  coor- 
dinates of  the  mid-point  of  the  S-curve.  Break  indicator  BR1  is  reset  to  the 
perimeter  length  from  the  starting  point  to  the  center  of  the  S-curve.  Vari- 
able DO  is  set  to  half  the  distance  from  the  starting  to  the  stopping  point. 

At  statement  50,  SGN  is  set  to  1 . If  the  shape  code  indicates  a 
left  circle  or  left  S-curve,  SGN  is  reset  to  -1.  Variable  V is  set  equal  to 
1-D/T0TLEN.  VS  is  the  square  of  V.  The  reciprocal  of  the  radius  of  curvature 
of  the  circle  or  the  circular  portion  of  the  S-curve  is  evaluated  using  a poly- 
nominal  approximation  to  the  solution  from  a transcendental  equation  containing 
the  reciprocal  of  the  radius  of  curvature.  The  approximate  radius  of  curva- 
ture, RPR,  is  improved  by  a series  of  linear  interpolations  if  the  value  for 
RPR  causes  an  error  greater  than  0.00001  in  the  transcendental  equation 

BR1*RPR  _ DD*RPR 
sin  2 Y~ 

When  RPR  is  within  the  desired  accuracy,  control  resumes  at  statement 
51.  The  radius  of  curvature,  R,  is  computed.  A temporary  variable,  ARG,  is 
evaluated.  The  height  of  the  center  of  the  circle  from  the  line  connecting  the 
starting  and  stopping  points,  H,  is  set  to  0.  If  variable  ARG  is  greater  than 
0,  H is  recomputed.  The  distance  to  the  first  break,  BR1 , is  tested  to  see 
whether  the  circular  arc  is  greater  than  half  a circle.  If  so,  the  sign  of 
the  height  is  changed.  The  x-  and  y-coordinates  of  the  center  of  the  circle 
are  computed.  The  components  of  the  vector  from  the  center  to  the  starting 
point.  Cl  1 and  C12,  are  computed.  All  variables  needed  to  compute  points  on 
the  S-curve  or  circle  are  now  available,  so  control  returns  to  the  calling 
program. 


Processing  continues  at  statement  60  for  the  remaining  shape  codes. 

If  the  total  perimeter  is  greater  than  the  straight-line  distance  from  start  to 
stop,  control  transfers  to  statement  65.  If  not,  and  if  variable  MODE  is  equal 
to  0,  a copy  of  the  input  card  is  printed.  An  error  message  is  also  printed 
indicating  that  the  map  distance  from  start  to  stop  exceeds  the  total  segment 
length.  The  shape  code  is  set  to  0,  and  control  returns  to  the  calling  pro- 
gram. 
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At  statement  65,  the  shape  code  is  tested;  if  neither  a right  nor  a 
left  rectangle  is  indicated,  control  transfers  to  statement  80.  Otherwise,  for 
a rectangular  segment,  the  distance  from  the  start  to  the  first  bend,  BR1 , is 
computed.  If  this  distance  is  greater  than  0.05  of  the  total  length,  control 
transfers  to  statement  70.  Otherwise,  the  rectangle  is  assumed  to  be  so  shal- 
low that  a straight-line  approximation  is  adequate,  and  the  shape  code  is  set 
to  0.  Control  then  returns  to  the  calling  program. 

At  statement  70  the  perimeter  to  the  second  bend  in  the  rectangle, 
BR2,  is  computed.  SX  and  SY,  the  x-  and  y-components  of  the  slope  of  the  vec- 
tor from  starting  point  to  stopping  point,  are  computed  and  control  returns  to 
the  calling  program. 

The  only  segments  that  reach  statement  80  are  the  angles.  The  first 
time  an  angle  shape  code  is  processed,  it  is  in  character  form  with  at  least 
the  first  two  character  positions  filled  with  binary  zeros.  After  this  shape 
code  has  been  processed  for  the  first  time,  it  will  be  replaced  by  a floating- 
point number  that  does  not  have  zeros  for  all  of  the  first  12  bits.  At  state- 
ment 80  control  transfers  to  statement  82  if  the  first  12  bits  of  the  shape 
code  are  zero.  If  not,  the  sign  of  the  shape  code  is  stored  in  variable 
SGN,  and  the  distance  to  the  vertex  of  the  angle  is  retrieved  as  the  magnitude 
of  the  shape  code  and  is  stored  in  variable  3R1.  Control  transfers  to  state- 
ment 140. 

At  statement  82,  SGN  is  set  to  0.  Variable  N is  set  to  0;  it  will 
contain  the  numerical  part  of  the  shape  code.  Variables  P 1 0 and  DPF  are  set 
to  1 . 

The  loop  through  statement  100  scans  the  characters  in  the  shape 
code.  The  characters  are  retrieved  in  variable  KAR.  If  no  character  is  found, 
control  goes  to  the  end  of  the  loop  and  another  character  is  sought.  If  the 
character  is  not  a decimal  point,  control  goes  to  statement  85.  If  a decimal 
point  is  found,  variable  DPF  is  set  to  10,  and  control  transfers  to  the  end  of 
the  loop. 


At  statement  85  the  character  is  checked  for  a digit.  If  no  digit 
is  found,  control  transfers  to  statement  90.  Otherwise,  the  cumulative  value 
of  the  string  is  stored  in  variable  N.  Variable  P10  is  multiplied  by  DPF , 
which  is  1 if  no  decimal  point  has  been  encountered,  or  10  if  a decimal  point 
has  been  encountered.  Control  transfers  to  the  end  of  the  loop. 

At  statement  90  a check  for  an  illegal  character  is  made.  If  an  il- 
legal character  is  found,  control  transfers  to  statement  110.  If  not,  SGN  is 
set  to  1.  If  the  character  is  an  L,  SGN  is  reset  to  -1.  If  the  loop  is  exited 
normally  after  the  shape  code  has  been  completely  scanned,  control  transfers  to 
statement  130. 

At  statement  110  the  card  image  is  printed  if  mode  is  equal  to  0.  An 
error  message  is  also  printed  indicating  the  illegal  character  in  th<.  shape 
code.  The  shape  code  is  set  to  0,  and  control  returns  to  the  calling  program. 

At  statement  130  the  distance  from  start  to  break,  which  is  part  of 
the  shape  code,  is  computed  and  converted  to  miles.  The  length  of  the  second 
leg  of  the  angle  is  computed  and  saved  in  variable  BR2.  A validity  checx  is 
performed  on  the  triangle  formed  by  the  two  sides  of  the  angle  and  the  line 
connecting  the  end  points.  If  any  sides  are  invalid,  that  is,  if  the  sum  of 
the  lengths  of  any  two  sides  is  less  than  the  length  of  the  third,  an  error 
message  is  printed.  Otherwise  control  transfers  to  statement  160.  If  an 
error  is  found,  the  shape  code  is  set  to  0 and  control  transfers  to  state- 
ment 20. 


At  statement  160  a temporary  variable,  F,  is  computed.  This  vari- 
able and  the  total  perimeter  of  the  angle  are  used  to  compute  the  height  of 
the  vertex  above  the  line  connecting  the  starting  and  stopping  nodes.  The 
height,  H,  is  then  used  in  the  computation  of  the  x-  and  y-coordinates  of 
the  vertex  of  the  angle.  Control  returns  to  the  calling  program. 

h.  Subroutine  COORD 

Subroutine  COORD  is  given  a distance,  in  miles,  from  the  beginning 
of  a segment  and  returns  the  coordinates  in  MCU.  Parameters  describing  the 
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segment  to  be  processed  have  been  stored  in  COMMON  block  COPARM  by  subroutine 
SHAPCOM  before  COORD  is  called.  Argument  CUMLEN  is  the  cumulative  length  along 
the  string,  in  miles;  arguments  XX  and  YY  are  the  coordinates  returned  for  a 
point  CUMLEN  miles  from  the  start  of  the  segment. 

The  first  statement  of  COORD  sets  S equal  to  the  cumulative  length. 

If  the  shape  code  is  nonzero,  control  transfers  to  statement  10.  The  coordi- 
nates of  the  point  on  a straight-line  segment  are  computed  and  returned  in 
variables  XX  and  YY.  Control  returns  to  the  calling  program. 

I 

At  statement  10  control  transfers  to  statement  30  if  the  shape  code 
indicates  other  than  a circular  or  S-curve  segment.  For  circular  and  S-curve 
segments,  the  reciprocal  of  the  radius  of  curvature  is  stored  in  RIP.  The  co- 
ordinates of  the  center  of  the  circular  portion  are  stored  in  XC  and  YC.  The 
components  of  the  vector  from  the  center  of  the  circle  to  the  initial  node  are 
stored  in  Cl  and  C2.  If  the  point  on  the  segment  is  less  than  or  equal  to 
0.999  of  the  first  break  distance  or  if  the  shape  code  indicates  a circular 
segment,  control  transfers  to  statement  20.  The  statements  following  this  test 
change  parameters  to  generate  coordinates  for  the  second  circular  portion  of  an 
S-curve.  The  sign  of  the  reciprocal  of  the  radius  of  curvature  is  reversed. 

The  cumulative  distance,  S,  is  set  to  the  distance  from  the  mid-point  of  the 
S-curve.  The  coordinates  of  the  center  of  the  second  circular  portion,  XC  and 
YC,  are  computed.  Variables  Cl  and  C2  are  recomputed  for  the  new  center. 

At  statement  20,  the  sine  and  cosine  of  the  angle  subtended  by  the 
perimeter  corresponding  to  S are  computed.  The  coordinates  XX  and  YY  of  the 
point  are  computed,  and  control  returns  to  the  calling  program. 

At  statement  30,  control  transfers  to  statement  60  if  the  shape  code 
indicates  that  the  segment  is  not  a rectangle.  Otherwise,  variable  SGN  is  set 
to  1 . If  the  shape  code  indicates  a left  rectangle,  SGN  is  reset  to  -1.  If 
S,  the  distance  along  the  rectangle,  is  greater  than  1.05  times  the  first 
side's  length,  control  transfers  to  statement  40.  If  S is  greater  than  0.95 
times  the  length  of  the  first  leg,  S is  set  to  the  length  of  the  first  leg. 

The  x-  and  y-coordi nates  of  the  point  on  the  first  leg  are  computed  by  linear 
interpolation,  and  control  returns  to  the  calling  program. 


33 


At  statement  40,  S is  tested  to  see  whether  it  falls  on  the  second 
leg  of  the  rectangle.  If  S is  greater  than  1.05  times  BR2,  the  length  of  the 
second  leg  of  the  rectangle,  control  transfers  to  statement  50.  If  S is 
greater  than  0.95  times  BR2,  S is  set  equal  to  BR2.  The  x-  and  y-coordinates 
of  the  point  on  the  second  leg  are  computed  by  linear  interpolation,  and  con- 
trol returns  to  the  calling  program. 

At  statement  50,  the  x-  and  y-coordinates  of  a point  on  the  third 
leg  of  the  rectangle  are  computed  by  linear  interpolation.  Control  returns  to 
the  calling  program. 

At  statement  60  the  distance,  S,  is  compared  to  the  length  of  the 
first  side  of  an  angle  segment.  IF  S is  greater  than  this  length,  control 
transfers  to  statement  70.  If  not,  the  x-  and  y-coordinates  are  computed  by 
interpolation  for  a point  on  the  first  leg.  Control  returns  to  the  calling 
program. 


At  statement  70  the  distance  along  the  angle  is  decreased  by  the 
length  of  the  first  leg  of  the  angle.  The  coordinates  of  the  point  on  the 
second  leg  are  computed  by  linear  interpolation,  and  control  returns  to  the 
calling  program. 

i . Subroutine  MAPPLT 

Subroutine  MAPPLT  draws  a map  of  the  street  segments,  one  line  per 
segment,  with  the  end-point  nodes  and  the  street  segments  numbered.  Up  to  10 
maps  can  be  drawn.  When  the  first  map  is  being  plotted,  MAPPLT  draws  all  of 
the  segments  on  a map-description  card  on  each  call.  After  the  first  map, 
MAPPLT  draws  the  entire  region  requested  on  one  call. 

Subroutine  MAPPLT  has  three  arguments.  Argument  II  indicates  the  se- 
quence number  of  the  map.  Arguments  KI  and  KF  are  the  numbers  of  the  initial 
and  final  segments  to  be  drawn  on  the  map  if  the  segments  lie  within  bounds. 

The  coordinates  of  the  region  bounding  the  map  are  contained  in  ar- 
rays in  COMMON  block  MPDATA.  In  this  COMMON  block,  arrays  XMIN  and  XMAX  are 
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the  minimum  and  maximum  x-coordi nates  for  the  map.  XLEN  is  the  length,  in 
inches,  of  the  map  in  the  x-direction.  YM1N , YMAX , and  YLEN  are  the  corre- 
sponding arrays  in  the  y-direction.  Array  YHCUT  contains  the  height,  in 
plotter  inches,  at  which  the  map  must  be  sliced  into  strips.  Array  SVAV  con- 
tains the  miles  per  MCU  conversion  factor  for  each  map.  Arrays  TRX  and  TRY 
contain  the  x-  and  y-components  of  the  translations  of  the  coordinate  systems 
of  the  maps  with  respect  to  the  coordinate  system  of  the  first  map.  Each  map 
has  its  own  coordinate  system.  Array  MSEQ  indicates  which  map  coordinate  sys- 
tem the  arrays  XMIN,  XMAX,  YMIN,  and  YMAX  are  in.  Variable  PLEN  is  the  length 
of  the  plot  in  plotter  inches.  It  is  the  total  length  from  the  start  of  the 
first  piece  to  the  end  of  the  last  strip  of  the  map.  Argument  CNVLEN  is  the 
miles  per  map  length  unit  conversion  factor. 

The  first  two  executable  statements  of  subroutine  MAPPLT  save  the 
initial  and  final  segment  numbers  in  variables  K1  and  K2.  If  the  map  number, 

II,  equals  1,  control  transfers  to  statement  20.  Otherwise,  the  ratio  of  miles 
per  MCU  conversion  factor  for  the  current  and  overall  maps  is  computed  and 
stored  in  variable  SCR.  The  x-  and  y-components  of  the  translation  of  the  cur- 
rent coordinate  system  relative  to  the  overall  coordinate  system  are  saved  in 
variables  TX  and  TY.  Since  all  segments  are  examined  after  the  first  map  has 
been  plotted,  variable  K1  is  set  to  1. 


At  statement  20,  if  the  first  map  is  being  drawn  and  subroutine 
MAPPLT  has  been  called  previously,  control  transfers  to  statement  110.  Other- 
wise, the  first-time-through  indicator  is  set  to  FALSE.  At  this  point  all 
parameters  that  will  apply  to  the  entire  map  are  set.  These  parameters  in- 
clude the  miles  per  MCU  conversion  (AVMD);  the  map  bounds  in  the  current  co- 
ordinate system  and  in  the  overall  coordinate  system;  the  height  of  a strip 
of  the  map;  the  maximum  length;  the  number  of  map  strips  (MX);  the  map  scale 
factors;  and  the  intervals,  in  MCU,  at  which  the  strips  are  cut.  These  pa- 
rameters are  printed  according  to  format  90. 

At  statement  100  the  segment  data  are  read  from  unit  1.  Unit  1 is 
tested  for  an  end-of-file.  No  end-of-file  will  be  encountered  unless  an  error 
exists  in  the  program.  When  no  end-of-file  is  encountered,  control  transfers 
to  statement  110. 
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At  statement  l'O  a loop  through  statement  200  tests  each  segment  to 
see  whether  it  falls  within  the  frame  of  the  map;  if  it  does,  it  will  be  plot- 
ted. Variables  NI  and  NF  are  set  equal  to  the  numbers  of  the  nodes  bounding 
the  segment.  The  mid-point  coordinates  of  the  segment  are  saved  in  variables 
XMD  and  YMD.  The  lines  in  the  node  number  array  at  which  the  initial  and 
final  nodes  occur  are  saved  in  variables  NS1  and  NS2.  The  coordinates 
of  these  nodes  are  retrieved. 

Initially  the  segment  is  assumed  to  be  entirely  within  the  bounds, 
and  indicators  I NB I , INBM,  and  INBF  are  set  to  1.  If  the  coordinates  of  the 
initial  node  lie  outside  the  frame  of  the  map,  INBI  is  set  to  0 Similar 
tests  are  made  on  the  coordinates  of  the  mid-point  of  the  segment  and  the  co- 
ordinates of  the  final  node  of  the  segment.  If  all  three  points  are  outside 
the  frame  of  the  map,  control  transfers  to  statement  200  and  the  segment  is 
not  plotted.  For  segments  that  are  at  least  partially  within  the  frame  of 
the  map,  the  shape  code  is  saved  in  variable  ISF.  If  the  shape  code  is  77  , 
the  segment  is  not  to  be  plotted  and  control  transfers  to  statement  200. 
Otherwise,  the  street  number  of  the  segment  and  its  total  length  in  miles  are 
saved  in  variables  NUMST  and  T0TLEN.  The  number  of  points  to  be  used  in  plot- 
ting half  the  segment  (NPMID)  is  computed.  The  number  will  be  restricted  to 
a maximum  of  10  points.  The  total  number  of  points  per  segment,  NPPSEG,  is 
set  to  twice  NPMID. 

Subroutine  SHAPC0M  is  called  to  set  up  the  parameters  needed  to  gen- 
erate coordinates  of  points  on  the  segment.  The  cumulative  length  along  the 
segment  is  initially  set  to  0.  A step  size,  DS,  is  computed  as  the  total 
length  divided  by  the  number  of  points  to  be  plotted  on  the  segment.  The 
coordinates  of  the  initial  node  are  converted  from  the  overall  coordinate 
system  to  the  current  coordinate  system  and  are  stored  in  variables  XX  and  YY. 
The  number  of  the  strip  of  the  map  into  which  the  node  falls  is  computed.  Both 
a current  value  of  the  strip  number,  NMAP,  and  a value  for  the  previous  point, 
NMAP0,  will  be  used.  The  pen  position,  up  or  down,  is  determined  by  whether 
the  initial  point  was  in  bounds.  Variable  IPEN  will  be  3 if  the  point  is  out 
of  bounds  and  2 if  the  point  is  in  bounds.  If  the  point  is  out  of  bounds,  con- 
trol transfers  to  statement  130.  If  not,  the  coordinates  of  the  point  are  con- 
verted to  plotter  inches  and  stored  in  variables  XP  and  YP.  If  the  current 
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node  has  already  been  plotted  as  the  last  node  on  the  previous  segment,  control 
transfers  to  statement  120.  If  not,  the  node  number  and  a small  square  marking 
its  position  are  appended  to  the  map.  At  statement  120  the  pen  is  moved  to  the 
position  of  the  current  point  on  the  segment. 

I 

Statement  130  starts  a loop  through  statement  170  that  will  advance 
the  pen  through  the  remaining  points  on  the  segment.  The  cumulative  length  is 
incremented  by  DS.  Subroutine  COORD  is  called  to  obtain  the  coordinates  of 
the  point  in  MCU. 

At  statement  140  the  coordinates  are  converted  to  plotter  inches. 

The  point  is  assumed  to  be  in  bounds,  and  variable  INB  is  set  to  1 . If  the 
coordinates  of  the  point  are  out  of  bounds,  INB  is  reset  to  0.  If  the  pen 
has  been  up  and  the  current  point  is  out  of  bounds,  or  if  the  strip  number  is 
greater  than  the  number  of  the  final  strip,  control  transfers  to  statement  160. 

Otherwise,  the  pen  is  moved  to  the  position  of  the  current  point.  If  the  pen 
is  up,  it  is  lowered.  Variable  IPEN  is  recomputed  to  reflect  whether  the  point 
is  in  bounds. 

At  statement  150,  if  the  loop  index  is  not  equal  to  the  number  of  the 
mid-point  of  the  segment,  control  transfers  to  statement  160.  Otherwise,  the 
segment  number  is  appended  to  the  map  near  the  segment  mid-point,  and  the  pen 
is  repositioned  at  the  mid-point. 

At  statement  160  the  number  of  the  current  strip  is  computed.  If  the 
current  strip  number  is  equal  to  the  previous  strip  number,  control  transfers 
to  statement  170.  If  not,  the  old  strip  number  (NMAP0)  is  set  equal  to  the 
current  strip  number;  IPEN  is  set  to  3,  indicating  that  the  pen  is  up;  and 
control  transfers  to  statement  140.  In  this  case,  the  pen  is  positioned  at 
the  current  point  on  the  new  strip. 

Statement  170  is  the  end  of  the  loop  that  causes  the  segment  to  be 
drawn.  If  the  last  point  drawn  is  out  of  bounds,  control  transfers  to  state- 
ment 200.  Otherwise  the  node  number  and  a small  square  marking  the  node's 
position  are  appended  to  the  map.  The  pen  is  repositioned  at  the  last  node. 

The  number  of  the  node  is  saved  in  variable  LASTNN. 
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Statement  200  is  the  end  of  the  loop  that  draws  the  various  seg- 
ments. At  statement  300,  if  the  map  being  drawn  is  other  than  the  first  map, 
the  plotter  pen  is  positioned  2 inches  beyond  the  end  of  the  last  strip. 
Control  returns  to  the  calling  program.  The  pen  is  moved  beyond  the  end  of 
the  first  map  in  the  main  program,  RCINPT. 

j.  Subroutine  LINEIN 

Subroutine  LINEIN  reads  information  from  card  images  in  free  for- 
mat. The  information  can  be  integer,  floating-point,  or  alphanumeric,  and 
alphanumerical  data  can  be  delimited  by  blanks,  asterisks,  dollar  signs,  or 
single  quotes. 

Subroutine  LINEIN  has  six  arguments.  Argument  I UN  is  the  number 
of  the  unit  from  which  the  card  image  is  to  be  read.  Argument  NIN  is  the 
number  of  data  items  to  be  read  from  the  card.  Argument  INPT  is  the  array 
into  which  the  data  are  placed. 

Argument  ITYPE  indicates  the  type  of  data  to  be  transferred.  Each 
octal  digit  of  ITYPE,  from  right  to  left,  specifies  the  type  for  the  data. 

If  the  digit  is  0,  an  integer  is  returned  by  LINEIN.  If  the  digit  is  1 , a 
floating-point  number  is  returned.  If  the  digit  is  2,  a word  of  a character 
string  delimited  by  *,  $,  or  1 is  returned.  If  the  digit  is  a 3 or  a 7,  a 
word  of  a blank-bounded  character  string  is  returned.  This  string  is  also 
terminated  by  an  equal  sign,  a left  parenthesis,  or  a slash.  Character 
strings  of  types  2 and  3 are  right-justified  with  preceding  binary  zeros; 
character  strings  of  type  7 are  left-justified  with  blank  fill. 

The  fifth  argument,  MODE,  causes  LINEIN  either  to  start  at  a new 
card  or  to  continue  reading  the  last  card.  A negative  MODE  causes  LINEIN  to 
continue  reading  the  last  card.  If  MODE  is  0,  reading  starts  at  a new  card. 
If  MODE  is  positive,  a new  card  is  read  and  printed. 

The  sixth  argument,  IBRK,  is  a break  and  error  indicator.  I5RK  is 
set  to  the  break  character  when  the  read  is  okay.  IBRK  is  set  to  0 when  a 
read  is  resumed  after  column  80,  to  -1  when  an  end-of-file  is  read,  or  to  -2 
when  an  error  is  detected  while  information  from  the  card  is  being  processed. 
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The  first  executable  statement  sets  logical  variable  DONE  to  FALSE. 
The  word  counter  (IW)  is  set  to  1 . If  MODE  is  less  than  0,  control  transfers 
to  statement  40.  If  not,  variable  IPRINT  is  set  to  MODE. 

At  statement  10  the  column  indicator  II  is  set  to  1.  The  break  in- 
dicator IBRK  is  set  to  -1.  A card  is  read  from  unit  IUN,  and  each  column  of 
the  card  is  placed  in  a separate  word  in  array  IC.  The  unit  is  tested  to  see 
whether  an  end-of-file  has  been  encountered  by  the  read.  If  so,  control  trans- 
fers to  statement  240.  If  not,  control  resumes  at  statement  30.  At  state- 
ment 30,  if  IPRINT  is  greater  than  0,  the  card  that  was  read  is  printed. 

At  statement  40,  a loop  through  statement  60  zeros  out  the  INPT  ar- 
ray for  the  variables  to  be  returned.  The  break  indicator  is  set  to  0.  If 
the  column  indicator,  II,  is  greater  than  80,  control  returns  to  the  calling 
program.  Otherwise,  variable  IAB  is  set  to  0.  IAB  is  used  to  hold  the  char- 
acter bounding  a bounded  text  string. 

The  loop  through  statement  65  presets  the  sign  and  magnitude  arrays 
ISGN  and  IV  to  be  used  in  building  numbers.  The  three  elements  in  each  array 
are  used  to  hold  the  integer  part,  the  fraction  part,  and  the  exponent  part  of 
a number.  The  sign  array  will  contain  either  +1  or  -1,  depending  on  the  sign 
of  the  corresponding  part  of  the  number.  The  magnitudes  of  the  three  parts 
are  initially  set  to  0.  Variable  ITYP  will  contain  the  digit  indicating  the 
type  of  the  word  being  processed.  Variables  IP  and  LB  are  set  to  1 . IP  indi- 
cates which  of  the  three  parts  of  the  number  is  being  processed.  LB  = 1 when 
the  previous  character  is  a blank.  Variable  NT  is  set  to  0.  It  will  hold  a 
count  of  the  number  of  digits  following  a decimal  point. 

The  loop  through  statement  200  controls  the  scan  of  the  card,  column 
by  column.  The  scan  begins  at  the  column  indicated  by  II  and  continues  through 
what  would  appear  to  be  column  81.  IC ( 81 ) is  preset  to  a blank  to  provide  a 
terminator  for  an  item  extending  into  column  80.  Variable  II  is  set  to  the 
number  of  the  next  column.  Variable  ICHAR  holds  the  character  in  the  current 
column.  If  ICHAR  is  a plus  sign  and  column  80  is  being  processed,  control 
transfers  to  statement  200,  the  end  of  the  loop.  If  not  and  if  ICHAR  is  a 
blank  (which  is  stored  in  variable  IBK),  and  the  previous  character  was  blank. 


and  no  bounded  character  string  is  in  process,  control  transfers  to  statement 
200.  If  not  and  if  IAB  is  0,  which  means  no  bounded  character  string  is  being 
processed,  control  transfers  to  statement  90.  Otherwise,  if  the  current  char- 
acter is  equal  to  the  break  character  on  a bounded  character  string,  or  if  the 
break  is  a blank  and  the  current  character  is  a comma  or  a right  parenthesis, 
control  transfers  to  statement  70.  Otherwise,  if  the  break  indicator  is  a 
blank  and  the  character  in  progress  is  a slash,  a left  parenthesis,  or  an  equal 
sign,  control  transfers  to  statement  68.  If  not,  the  present  character  is  ap- 
pended to  the  word  of  text  being  built. 

The  word  of  text  is  kept  in  variable  II.  The  present  character  is 
added  to  the  end  of  this  word,  and  counter  IAN  is  incremented  by  1 . If  IAN 
is  less  than  or  equal  to  10,  control  transfers  to  statement  200.  Otherwise, 

IAN  is  set  to  1 and  control  transfers  to  statement  159,  where  the  text  in  II 
will  be  moved  to  the  INPT  array. 

At  statement  68,  if  IAN  is  equal  to  1,  control  transfers  to  state- 
ment 220.  Otherwise,  the  word  of  text  being  built  in  II  must  be  moved  to  the 
INPT  array.  In  this  case,  variable  DONE  is  set  to  TRUE. 

At  statement  70  the  break  indicator  is  set  to  0.  If  a new  word  of 
text  has  not  been  started,  control  transfers  to  statement  200.  Otherwise,  if 
the  type  of  the  text  is  not  type  7,  control  transfers  to  statement  159.  Other- 
wise, blanks  are  appended  to  the  end  of  the  word  being  built  in  II,  and  con- 
trol transfers  to  statement  159. 

At  statement  90,  LB  is  set  to  0 because  the  current  character  is 
not  blank.  If  the  character  is  not  a digit,  control  transfers  to  statement 
100.  If  the  character  is  a digit,  it  is  changed  to  its  binary  value,  and  the 
cumulative  value  of  the  number  being  processed  is  stored  in  variable  IV.  If 
a decimal  point  has  been  encountered,  variable  NT  is  incremented  by  1.  Con- 
trol transfers  to  statement  200. 

At  statement  100  a test  is  made  for  an  exponent  field.  If  the  pres- 
ent character  is  not  an  E,  or  if  the  type  of  the  number  is  not  integer  or 
floating-point,  or  if  an  exponent  field  has  already  been  processed,  control 
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transfers  to  statement  102.  Otherwise,  IP  is  set  to  3,  and  control  transfers 
to  statement  200. 

At.  statement  102,  a loop  through  statement  130  is  executed.  In  this 
loop,  the  present  character  is  tested  for  one  of  12  special  characters.  If 
the  character  is  not  one  of  these  12,  control  transfers  to  the  end  of  the 
loop.  If  one  of  the  12  is  found,  control  transfers  to  a statement  that  does 
the  appropriate  processing.  If  a dollar  sign,  an  asterisk,  or  a single  quote 
is  found,  indicating  the  start  of  a character  string,  control  transfers  to 
statement  105.  Otherwise,  if  the  type  of  the  word  is  not  equal  to  2,  control 
transfers  to  statement  131,  where  an  error  message  is  printed.  When  the  type 
is  2,  IAN  is  set  to  1,  indicating  that  the  next  character  position  to  be  filled 
is  the  first  character.  IAB  is  set  equal  to  the  present  character,  which  will 
terminate  the  string  at  its  next  occurrence.  Control  transfers  to  statement 
200.  If  a decimal  point  is  found,  control  transfers  to  statement  110. 

At  statement  110  variable  IP  is  incremented  by  1.  This  causes  the 
processing  to  build  I V ( 2 ) as  the  value  of  the  number  following  the  decimal 
point.  If  IP  > 2,  control  transfers  to  statement  131  where  an  error  message 
is  printed.  Otherwise,  control  transfers  to  statement  200. 

At  statement  120  a minus  sign  is  processed.  If  the  sign  of  the  num- 
ber being  built  is  already  negative,  or  if  the  fractional  part  of  a number  is 
being  built,  control  transfers  to  statement  131.  Otherwise,  the  sign  of  this 
part  of  the  number  is  set  to  -1,  and  control  transfers  to  statement  200. 

At  statement  125  the  left  parenthesis,  slash,  and  equal  sign  charac- 
ters are  processed.  Variable  DONE  is  set  to  TRUE.  If  the  type  of  the  word  is 
0 or  1 , control  transfers  to  statement  140,  where  a numeric  value  will  be  pro- 
cessed. Otherwise,  control  transfers  to  statement  159,  where  a word  of  text 
will  be  processed. 

Statement  130  is  the  end  of  the  loop  that  searches  for  the  12  special 
characters.  Any  character  that  is  not  processed  by  the  end  of  the  130  loop  is 
assumed  to  be  the  first  character  of  a string  of  blank-bounded  text.  If  the 
type  is  not  3 or  7,  or  if  the  character  is  not  alphanumeric,  control  transfers 
to  the  error  message  at  statement  131.  When  valid  alphanumeric  characters  are 


encountered,  IAN  is  set  to  2,  indicating  that  the  next  character  to  be  pro- 
cessed will  be  stored  in  the  second-character  position.  IAB  is  set  to  a 
blank.  The  present  character  is  stored  II,  which  will  eventually  hold  the 
word  of  text.  Control  transfers  to  statement  200. 

At  statement  131  an  error  message  is  printed  indicating  which 
column  of  the  card  image  contains  the  error.  The  card  image  is  also  printed. 
The  error  indicator,  IBRK,  is  set  to  -2,  and  control  returns  to  the  calling 
program. 


At  statement  140,  the  final  evaluation  of  numeric  constants  is 
performed.  F2,  the  fractional  part  of  the  number,  is  initially  set  to  0. 

If  NT  is  greater  than  0,  F2  is  re-evaluated.  If  the  type  does  not  indicate 
a floating-point  number,  control  transfers  to  statement  150.  Otherwise,  FI 
is  evaluated  as  the  floating-point  value  of  the  integer  that  is  to  be  re- 
turned. The  integer  portion  of  FI  is  stored  in  II,  and  control  transfers  to 
statement  159. 

At  statement  150,  the  desired  floating-point  number  is  evaluated 
and  stored  in  FI . 

At  statement  159,  the  value  in  II  is  transferred  to  the  INPT 
array.  Note  that  variables  FI  and  II  are  equi valenced.  The  word  counter, 

IW,  is  incremented  by  1 . If  IW  is  greater  than  the  number  of  words  to  be 
returned  by  LINE  IN , or  if  logical  variable  DONE  is  TRUE,  control  transfers 
to  statement  220.  Otherwise,  the  type  of  the  next  word  to  be  returned  is 
stored  in  ITYP. 

The  loop  through  statement  170  resets  the  sign  and  magnitude  arrays 
used  to  hold  the  appropriate  parts  of  numeric  constants.  II  and  NT  are  reset 
to  0.  LB  and  IP  are  reset  to  1.  The  loop  that  processes  the  columns  of  the 
card  terminates  at  statement  200.  If  column  80  of  the  card  is  a plus  sign, 
control  transfers  to  statement  10,  where  another  card  will  be  read.  Other- 
wise, if  the  next  column  to  be  processed  is  equal  to  the  character  bounding 
a text  string,  the  column  counter  II  is  incremented  by  1.  The  break  charac- 
ter is  returned  in  variable  IBRK. 
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At  statement  240  control  returns  to  the  calling  program. 


An  entry  point,  PRNTC,  is  provided  to  print  the  current  card  image. 
The  card  image  is  printed,  and  control  returns  to  the  calling  program. 

k.  Program  RCINPT 

Main  program  RCINPT  controls  the  reading  of  the  map  data  and  the 
plotting  of  the  maps.  It  uses  six  files:  INPUT,  OUTPUT,  TAPE 1 , TAPE2,  TAPE3 , 
and  TAPE8.  File  TAPE5  is  equivalenced  to  INPUT  in  order  to  test  for  end-of- 
records.  File  TAPE8  is  used  for  Calcomp  plot  output.  Files  TAPE1 , TAPE2,  and 
TAPE3  are  used  for  segment,  node,  and  street-name  data,  respectively. 

Blank  COMMON  contains  a title  array  and  storage  for  the  segment  data. 
COMMON  block  COPARM  has  been  described  under  subroutine  SHAPCOM.  COMMON  block 
NDDATA  contains  the  node  data.  The  items  in  this  COMMON  block  are  a count  of 
the  number  of  nodes  (KNODES),  an  array  of  the  numbers  of  the  segments  bounding 
the  node  (NBS),  the  node  number  (NODNUM),  the  number  of  times  the  node  has  oc- 
curred (TIMNOD),  and  the  x-  and  y-coordinates  of  the  node  (XNOD  and  YNOD). 
COMMON  block  MPDATA  contains  map  information  and  is  described  under  subroutine 
MAPPLT.  Array  INPT  holds  data  returned  by  subroutine  LINEIN.  Since  these  can 
be  either  integer  or  floating-point  data,  array  FNPT  is  equivalenced  to  INPT 
so  that  the  different  types  of  data  can  be  accessed  properly.  Arrays  ISTG  and 
STG  are  equivalenced  primarily  for  convenience  in  storing  and  retrieving  inte- 
ger and  floating-point  segment  data.  Variable  XNI  is  equivalenced  to  XN ( 1 ) and 
variable  YNI  to  YN(1)  primarily  for  convenience  in  writing  certain  FORTRAN 
sta tements. 

Several  variables  are  preset  in  DATA  statements.  Variable  CMDMXR, 
the  maximum  error  in  map  distance  conversion,  is  preset  to  .1.  When  the  rela- 
tive error  in  map  distance  conversion  exceeds  this  number,  the  program  causes 
the  map  distance  conversion  to  be  redefined  in  one  mode  of  operation  or  prints 
a warning  message  in  the  other.  Variable  MAXSEG  gives  the  maximum  number  of 
segments  that  may  be  used.  It  must  be  equal  to  the  second  dimension  of  array 
STG.  For  clarity,  the  first  subscript  of  the  segment  data  array,  STG,  is  used 
in  symbolic  form  so  that  the  particular  type  of  segment  data  will  be  apparent. 
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The  numerical  values  for  the  symbolic  names  are  given  in  a DATA  statement. 
Variable  NSTR  indicates  the  subscript  for  street  numbers.  Variables  NN 1 and  NN2 
indicate  the  initial  and  final  node  numbers.  Variable  LEN  is  the  subscript  for 
the  segment  length.  NH  is  the  subscript  for  the  number  of  houses  on  the  seg- 
ment. NSPD  is  the  subscript  for  the  speed  limit  on  the  segment.  NWAY  is  the 
subscript  for  the  number  of  ways  the  street  may  be  traveled.  NRQF  is  the  sub- 
script for  the  refuse-quantity  adjustment  factor.  NXMID  and  NYMID  are  the  sub- 
scripts for  the  x-  and  y-coordinates  of  the  segment  midpoint.  NSF  is  the  sub- 
script for  the  shape  code  of  the  segment.  Variable  MODE  controls  the  printing 
of  the  input  data  map-description  cards.  When  MODE  is  1,  the  data  cards  are 
printed  as  they  are  read  by  LINE  IN . The  total  number  of  houses  and  total  ref- 
use quantity  are  set  to  0 in  a DATA  statement.  Arrays  holding  the  map  distance 
conversion  factors  (SVAV)  and  the  x-  and  y-coordinates  of  map  translations  (TRX 
and  TRY)  are  also  set  to  0 in  a DATA  statement. 

The  first  executable  statement  initializes  the  plotting  package.  The 
pen  is  moved  down  and  then  up  3 inches  to  force  at  least  a 3-inch  border  at 
the  bottom  of  the  plot.  The  input  unit  number  (IUN)  is  set  to  5.  This  number 
will  be  used  in  calls  to  subroutine  LINEIN.  Variable  IQUIT  is  set  to  0.  It 
will  hold  a count  of  errors  found  in  the  map-description  cards.  The  variables 
that  hold  counts  of  the  segments  and  nodes,  KF  and  KNODES,  are  set  to  0.  The 
maximum  number  of  errors  allowed  in  map-description  data  cards  before  proces- 
sing terminates  is  set  to  20. 

Subroutine  STRINP  is  called  to  input  the  street-name  data.  Argument 
N 1 1 R is  1 if  street-name  data  are  absent,  or  2 if  street-name  data  are  present. 
In  normal  operation,  street-name  data  will  be  present.  The  double  loop  through 
statement  10  zeros  out  the  segment-data  array,  ISTG.  The  loop  through  state- 
ment 30  reads  up  to  10  map-bounds  description  cards.  The  loop  is  executed  11 
times,  at  most,  to  ensure  that  an  end-of-record  indicator  will  be  encountered 
after  a tenth  data  card. 

XMIN,  XMAX,  YMIN,  and  YMAX  are  the  minimum  and  maximum  x-  and  y- 
coordinates  bounding  the  map.  Variables  XLEN  and  YLEN  are  the  x-  and  y-lengths 
of  the  map  in  plotter  inches.  Variable  YHCUT  is  the  height,  in  inches,  at  which 
the  map  should  be  cut  into  strips.  YHCUT  will  be  set  equal  to  30  inches  at 
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statement  25  if  this  field  is  zero  or  blank  on  the  data  cards.  Variable  MSEQ 
is  the  sequence  number  of  the  coordinate  system  for  arguments  XMIN,  XMAX,  YMIN, 
and  YMAX . If  MSEQ  is  blank  or  zero  on  the  data  card,  it  will  be  set  to  1. 

The  loop  is  executed  one  time  more  than  the  number  of  map-data  cards 
so  that  the  end-of-record  card  will  be  passed.  At  statement  40  the  number  of 
maps  (MAPS)  is  set  to  1-1. 

If  any  maps  are  to  be  drawn,  subroutine  MAPGRID  is  called  to  draw 
the  appropriate  grid  for  the  first  map.  Scale  ratio,  SCR,  is  set  to  1.  Vari- 
able MDFILE,  the  number  of  the  current  data  file,  is  set  to  0.  The  statements 
through  this  point  are  executed  only  once. 

At  statement  50  the  reading  of  the  map-description  data  begins.  Con- 
trol will  return  to  statement  50  each  time  a map  has  been  completely  processed. 
The  first  two  cards  in  the  map  description,  which  contain  scale  factors  and  de- 
fault values,  are  read  according  to  format  60.  UNLEN  is  the  unit  of  map  length 
measurement.  SCALEM  is  the  map  scale  in  feet  per  inch.  SCALEC  is  the  coordin- 
ate scale  in  plotter  inches  per  MCU.  CTYPE  is  the  cocrdi nate-use  type.  The 
default  speed  limit,  number  of  ways  of  traveling  a street,  number  of  sides  col- 
lected on  one  pass,  and  refuse  quantity  adjustment  factor  (RQAF)  are  read  from 
the  second  card  and  are  saved  in  variables  SPDDEF,  NWAYDEF,  NSIDDEF,  and  RQFDEF. 
If  an  end-of-file  or  end-of-record  is  encountered,  control  transfers  to  state- 
ment 1000;  otherwise  execution  continues  at  statement  70. 

At  statement  70  the  count  of  map  data  files  is  incremented  by  1 . If 
the  coordinate  scale  is  less  than  or  equal  to  0,  it  is  changed  to  1 inch  per 
MCU.  AVMD  and  AVMDDEF  are  the  map  distance  conversions  in  miles  per  MCU. 

CCMD,  the  count  of  strings  used  in  the  evaluation  of  AVMD,  is  preset  to  0. 

If  AVMD  is  not  0,  CCMD  is  reset  to  1 . If  the  map  scale  is  less  than  or  equal 
to  0,  it  is  reset  to  400  feet  per  plotter  inch.  The  map  length  conversion 
is  set  to  0.01  miles  per  map  length  unit.  If  UNLEN  has  been  specified  as 
positive,  the  conversion  is  re-evaluated.  The  coordinate-use  type  (CTYPE) 
is  tested  and  will  end  by  being  equal  to  either  the  letters  AVG,  FIRST,  or 
LAST.  The  default  speed  limit  is  set  to  15  miles  per  hour  if  the  default  on 
the  card  is  not  positive.  The  default  number  of  ways  of  traveling  a street  is 
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set  to  2 if  no  positive  number  is  found  on  the  card.  The  default  number  of 
sides  collected  on  one  pass  is  set  to  2 if  a positive  number  is  not  found.  The 
default  RQAF  is  set  to  1 if  a positive  factor  is  not  found  on  the  card.  Logi- 
cal variable  FIRSTT  is  evaluated  and  will  have  the  value  TRUE  the  first  time  a 
card  is  read.  If  the  first  map-description  record  is  being  processed,  or  if 

AVMU  is  positive,  control  transfers  to  statement  90.  Otherwise,  an  error  mes- 

sage is  printed  indicating  that  only  one  map  is  allowed  when  the  variable  map- 
coordinate  option  is  used;  the  job  is  then  terminated. 

At  statement  90  the  map  distance  conversion  is  saved  in  array  SVAV. 

On  all  maps  after  the  first,  the  scale  ratio,  SCR,  is  evaluated.  A count  in- 
crement, KINC,  is  set  to  0.  The  title  of  the  problem  and  a number  of  vari- 

ables describing  the  map  are  printed.  This  information  is  useful  primarily 
when  the  output  map  differs  drastically  from  what  was  desired. 

At  statement  210,  the  reading  of  the  map-description  cards  is  begun. 
The  street  number  and  the  number  of  the  first  node  are  read.  IBRK  is  tested 
to  see  whether  the  card  contains  an  error.  If  so,  control  transfers  to  state- 
ment 220.  If  not,  IBRK  is  tested  for  three  conditions:  an  end-of-file,  a call 
to  LINEIN  with  the  card  column  pointer  beyond  column  80,  or  a valid  read.  If 
an  end-of-file  (or  end-of-record)  is  encountered,  control  transfers  to  state- 
ment 50,  where  the  first  card  in  a new  map-description  record  will  be  read. 

It  should  be  impossible  for  a card  to  be  started  past  column  80. 

At  statement  220,  if  MODE  is  0,  the  card  being  processed  is  printed. 
An  error  message  is  then  printed  indicating  that  the  card  contains  some  error. 
The  error  counter  (IQUIT)  is  incremented  by  1 . If  the  number  of  errors  does 
not  yet  exceed  the  maximum,  control  transfers  to  statement  210  and  another 
map-description  card  is  read.  Otherwise,  the  job  terminates. 

Control  resumes  at  statement  240  if  no  errors  have  been  encountered 
during  the  reading  of  the  street  number  and  first  node  number.  The  count  of 
segments  is  incremented  by  1.  The  number  of  the  street,  if  present,  is  stored 
in  NUMST ; otherwise,  NUMST  is  set  to  0.  The  initial  node  number  is  saved 
in  NODE  I and  is  also  stored  in  the  segment  data  table.  Variable  T0TLEN,  a 
total  length  measurement,  is  set  to  0.  Function  IFIND  is  used  to  find 
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the  initial  node  number  in  the  node  data  table.  If  the  node  is  present  in  the 
table,  NS  1 will  be  positive,  and  control  transfers  to  statement  245.  If  not, 
NS1  is  negative  and  indicates  where  in  the  node  table  the  node  number  should  be 
stored.  Its  sign  is  changed,  subsequent  entries  in  the  table  are  moved  down  to 
make  room  for  the  new  node  number,  and  the  node  number  is  stored  in  the  table. 

A count  of  the  number  of  nodes  is  incremented  by  1 . At  statement  245  the  line 
number  in  the  table  is  saved  in  variable  NS2. 

At  statement  250,  LINEIN  is  called  to  obtain  the  segment  length,  num- 
ber of  houses  on  the  right  side  of  the  street,  number  of  houses  on  the  left 
side  of  the  street,  and  final  node  number  of  the  segment.  The  break  indicator 
is  tested.  Control  transfers  to  statement  220  if  an  error  exists,  to  state- 
ment 270  if  the  read  resumed  past  column  80,  or  to  statement  260  if  the  read 
was  valid. 


At  statement  260  the  first  number  returned  by  LINEIN  is  tested.  If 
it  is  zero,  no  data  were  found  on  the  card  (zero  segments  lengths  are  not  al- 
lowed), and  control  transfers  to  statement  270.  When  a valid  segment  length 
is  found,  the  segment  number  is  saved  in  array  NBS  indicating  that  the  segment 
bounds  the  most  recently  encountered  node.  Other  information  about  the  seg- 
ment is  stored  in  the  segment  data  array  (STG),  and  the  total  length  of  the 
present  string  is  incremented  by  the  length  of  the  current  segment.  The  num- 
ber of  houses  on  the  left  side  of  the  street  is  saved  in  the  NHL  array.  The 
final  node  number  is  saved  in  the  segment  data  array  and  in  variable  N0DEF. 

The  segment  counter  (KF)  is  incremented  by  1 . If  the  final  node  number  is  not 
positive,  control  transfers  to  the  error  message  at  statement  220.  Otherwise, 
the  second  node  number  is  sought  in  the  node  number  table.  If  it  is  present, 
NS2  will  indicate  its  line  number,  and  control  transfers  to  statement  265.  If 
NS2  is  not  positive,  it  indicates  where  the  node  should  be  stored  and  the  node 
is  appended  to  the  node  data  table. 

At  statement  265  the  segment  number  is  saved  in  the  NBS  array,  indi- 
cating that  the  segment  bounds  the  current  node.  The  break  character  returned 
by  LINEIN  is  tested.  If  it  is  neither  a slash  nor  a left  parenthesis,  control 
transfers  to  statement  250  where  another  segment  will  be  read.  If  a slash  or 
a left  parenthesis  is  encountered,  the  segment  count  is  decremented  by  1 at 
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statement  270.  Because  the  addition  of  the  terminal  node  number  to  the  node 
number  table  may  have  changed  the  location  of  the  initial  node  number,  func- 
tion IFIND  is  used  to  again  find  the  location  of  the  first  node. 


The  first  four  locations  of  the  I NPT  array  are  set  to  0.  If  the 
break  character  of  the  last  call  to  LINEIN  was  not  a left  parenthesis,  then 
the  speed  limit,  one-way  indicator,  number  of  sides  collected  on  one  pass,  and 
RQAF  may  be  on  the  card  before  the  coordinates  of  the  first  node.  In  this  case, 
LINEIN  is  called  to  read  these  items.  If  an  error  is  encountered,  control 
transfers  to  statement  220. 

At  statement  280  the  four  items  are  stored  in  the  appropriate  vari- 
ables; if  any  of  the  four  items  is  0,  the  default  values  are  used.  The  loop 
ending  on  statement  290  appends  to  the  segment  data  the  speed  limit,  street 
number,  number  of  ways  of  travel,  number  of  houses  to  be  serviced,  and  RQAF.  A 

count  of  the  total  number  of  houses  and  the  total  amount  of  refuse  is  accumu- 

lated in  this  loop.  If  the  break  character  returned  by  LINEIN  is  a left  pa- 
renthesis, control  transfers  to  statement  310,  where  the  coordinates  of  the 
first  node  will  be  processed.  Otherwise,  the  shape  code  SF  is  set  to  0.  If 
the  break  indicator  shows  that  the  read  resumed  after  column  80,  or  if  the 
break  character  is  a blank,  control  transfers  to  statement  340.  Otherwise, 
an  error  message  is  printed  indicating  an  illegal  character  in  the  card  just 
read.  Control  transfers  to  statement  210,  where  the  next  card  is  read. 

At  statement  310,  LINEIN  is  called  to  read  the  coordinates  of  the 
first  node.  If  an  error  is  encountered,  control  transfers  to  statement  220. 

If  not,  variable  NRJMP  is  used  when  a new  map  is  started  to  indicate  whether 

the  starting  or  ending  node  matches  a node  on  a previous  map.  This  variable  is 

set  to  1 at  present.  The  shape  code  is  obtained  from  the  INPT  array.  If  the 
first  two  numbers  returned  by  LINEIN  indicate  no  coordinates  or  zero  coordin- 
ates for  the  first  node,  control  transfers  to  statement  319,  where  the  coor- 
dinates of  the  ending  node  will  be  processed.  Otherwise,  the  coordinates  of 
the  first  nodes  are  saved  in  variables  SVX  and  SVY.  If  the  first  map- 
description  card  of  a record  is  being  processed  and  the  present  node  has  not 
been  used  before,  control  transfers  to  statement  320.  If  the  present  card  is 
not  the  first  map-descri ption  card  in  a record,  control  transfers  to  statement 
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315.  Otherwise,  translations  for  this  map  relative  to  the  original  map  are 
computed.  The  translations  are  printed.  Logical  variable  FIRSTT  is  set  to 
FALSE. 


At  statement  315,  if  the  coordinate-use  type  is  LAST,  the  count  of 
occurrences  of  this  node  is  set  to  0.  If  averaged  coordinates  are  not  to  be 
used,  and  if  the  count  of  node  occurrences  is  greater  than  0,  control  trans- 
fers to  statement  319.  Otherwise,  the  count  of  occurrences  of  the  node,  TIMNOD, 
is  incremented  by  1;  and  the  coordinates,  either  initial  or  averaged,  are  saved 
in  arrays  XNOD  and  YNOD. 

At  statement  319,  NRJMP  is  set  to  2.  At  statement  320,  if  the  last 
break  character  is  not  a left  parenthesis,  control  transfers  to  statement  340. 
Otherwise,  LINEIN  is  called  to  obtain  the  coordinates  of  the  final  node.  If 
an  error  is  encountered,  control  transfers  to  statement  220.  If  the  read  re- 
sumed after  column  80,  no  coordinates  are  present  and  control  transfers  to 
statement  340.  If  the  read  was  valid,  control  resumes  at  statement  330.  If 
the  coordinates  are  zero,  control  transfers  to  statement  340.  Otherwise,  the 
succeeding  statements  through  339  process  the  coordinates  of  the  ending  node 
in  the  same  manner  as  that  used  to  process  the  coordinates  of  the  first  node. 

At  statement  339,  if  NRJMP  is  equal  to  1,  control  transfers  back  to 
statement  315.  At  statement  340,  if  the  card  being  processed  is  not  the  first 
map-descri ptior,  card  in  a record,  control  transfers  to  statement  342.  Other- 
wise, an  error  message  is  printed  indicating  that  the  string  does  not  start  or 
end  on  a previously  defined  node.  The  error  counter,  I QUIT , is  incremented, 
and  control  resumes  at  statement  342. 

A loop  tnrough  statement  370  verifies  that  coordinates  have  been 
given  for  the  initial  <nd  final  nodes  of  the  string.  If  coordinates  have  not 
been  given  for  one  or  the  other,  an  error  message  is  printed. 

Following  statement  380,  if  the  shape  code  is  nonzero,  control 
transfers  to  statement  430.  Otherwise,  the  string  of  street  segments  is 
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straight.  Since  only  straight-segment  strings  can  be  used  to  redefine  a map 
scale,  testing  begins  for  a new  map-scale  definition.  Variable  CONMD  is  eval- 
uated as  the  scale  in  miles  per  MCU  for  the  current  string.  If  CCMD  is  0, 
control  transfers  to  statement  410.  Otherwise,  the  relative  error  between  the 
old  and  current  map  scales  is  computed.  If  AVMDDEF  is  0,  the  variable  map- 
scale  option  is  in  use  and  control  transfers  to  statement  400.  If  the  relative 
error  is  less  than  that  caused  by  a two-unit  error  in  the  string,  control  trans- 
fers to  statement  430.  Otherwise,  an  error  message  is  printed  indicating  that 
the  previous  string  deviates  badly  from  the  default  scale.  Control  transfers 
to  statement  430. 

At  statement  400,  if  the  relative  error  is  less  than  a permissible 
amount,  control  transfers  to  statement  430.  (The  permissible  amount  is  about 
10  percent  fur  long  strings  and  increases  to  40  percent  for  a string  of  one 
length-unit.)  Otherwise,  CCMD,  the  count  of  strings  used  in  the  current  scale 
definition,  is  reset  to  0.  A message  is  printed  indicating  a change  in  map 
scale,  and  ‘-.he  new  scale  is  computed. 

At  statement  430,  subroutine  SHAPC0M  is  called  to  process  the  pres- 
ent shape  code. 

At  statement  600  variable  CUMLEN  is  set  to  0.  The  loop  through  state- 
ment 690  will  plot  the  various  segments  on  the  map-description  string  just  read. 
The  loop  index  (K)  is  both  the  segment  number  and  its  line  number  in  the  seg- 
ment table.  The  segment  length  is  retrieved  in  variable  SLEN.  Subroutine 
COORD  is  called  to  obtain  coordinates  of  the  midpoint  of  the  segment.  If  the 
shape  code  is  neither  circular  nor  straight,  control  transfers  to  statement  610. 
Otherwise,  the  shape  code  is  appended  to  the  segment  data  table.  Control  trans- 
fers to  statement  650. 

At  statement  610,  if  the  shape  code  indicates  a rectangle,  control 
transfers  to  statement  630.  If  the  shape  code  indicates  a right  or  left  S- 
curve,  control  transfers  to  620.  Otherwise,  the  remaining  angle  shape  codes 
are  processed  by  the  statements  following  statement.  610.  A zero  is  stored  in 
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the  segment  data  table  as  the  shape  code  for  any  straight  segments  that  do  not 
include  the  vertex  of  the  angle.  A floating-point  number  is  stored  as  the 
shape  code  for  the  segment  including  the  vertex.  This  number  indicates  the 
length  from  the  start  of  the  segment  to  the  vertex  of  the  angle.  Control 
transfers  to  statement  650. 

At  statement  620,  segments  on  an  S-curve  are  processed.  Circular 
portions  other  than  the  center  of  the  S-curve  have  a circular  shape  code  stored 
in  the  segment  table.  If  a segment  includes  the  center  of  the  S-curve,  an  S- 
curve  shape  code  is  stored  for  it  in  the  segment  data  table.  Control  transfers 
to  statement  650. 

At  statement  630  rectangular  strings  are  processed.  A zero  is  stored 
in  the  segment  data  table  as  the  shape  code  for  straight  segments.  If  either 
corner  of  the  rectangle  falls  on  the  interior  of  a segment,  an  angle  shape  code 
is  stored  for  the  segment.  If  both  corners  fall  on  one  segment,  a rectangle 
shape  code  is  stored. 

Following  statement  650  the  cumulative  length  along  the  string  is 
evaluated.  If  the  segment  is  the  last  of  the  string,  control  transfers  to  the 
end  of  the  loop  at  statement  690.  Otherwise,  the  coordinates  of  the  ending 
node  of  the  segment  are  found  by  subroutine  COORD;  if  they  have  not  yet  been 
defined  or  if  average  coordinates  are  desired,  they  are  stored  in  the  XN0D  and 
YN0D  arrays. 

Statement  690  terminates  the  loop  that  evaluates  shape  codes  for  each 
segment  in  the  map-description  string.  If  two  sides  of  the  street  are  collec- 
ted on  one  pass,  control  transfers  to  statement  800.  If  collection  is  from  one 
side  of  a street  at  a time,  each  segment  in  the  map-description  string  just  pro- 
cessed will  be  repeated  as  a one-way  segment  in  the  opposite  direction,  with 
the  number  of  houses  on  the  left  side  of  the  street  in  the  segment  data  table. 
The  number  of  segments  to  be  repeated  is  evaluated  in  variable  KINC. 

The  loop  through  statement  700  repeats  the  previous  KINC  segments, 
but  makes  them  one-way  in  the  opposite  direction.  The  shape  codes  are  set  to 

77  , indicating  that  no  additional  plotting  is  necessary  on  these  segments. 
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At  statement  800  subroutine  MAPPLT  is  called  to  plot  the  segments 
from  the  map-description  string  just  processed.  The  count  of  segments  is  in- 
cremented by  KINC,  and  KINC  is  reset  to  0.  Control  transfers  to  statement 
210  for  the  reading  of  another  map-description  string. 

Statement  1000  is  reached  when  an  end-of-record  is  encountered  as 
a map-description  record  is  started.  The  plotter  pen  is  advanced  beyond  the 
end  of  the  first  map.  The  problem  title  and  the  number  of  houses  and  total 
refuse  quantity  are  printed  on  the  output.  If  no  map-description  strings  have 
been  encountered,  control  transfers  to  statement  1500.  Otherwise,  a list  of 
the  segments  is  printed  using  format  1010.  The  title  is  reprinted  at  the  top 
of  another  page,  and  the  node  data  are  printed.  The  segment  data  are  written 
to  file  1,  and  the  node  data  are  written  to  file  2.  An  end-of-file  is  written 
on  each  file.  If  one  or  no  map  was  plotted,  control  transfers  to  statement 
1500. 


The  loop  through  statement  1420  plots  the  grid  and  map  for  each  mao 
after  the  first.  A message  is  printed  to  show  that  the  map  has  been  plotted. 

At  statement  1500  an  additional  end-of-file  is  placed  on  the  plot 
tape.  The  terminating  call  to  the  plot  package  is  then  executed.  System 
subroutine  EXIT  is  called,  and  control  returns  to  the  system. 


52 


SECTION  IV 
INPUT  AND  OUTPUT 


1.  INPUT 

The  only  input  to  program  RCINPT  is  card  input.  Since  the  preparation 
of  data  cards  for  program  RCINPT  is  described  in  a separate  report  (Refer- 
ence 1),  only  the  cards  will  be  described  here. 

In  this  report,  the  word  record  means  all  data  cards  between  two  end-of- 
record  indicators  (cards  with  7-8-9  multipunched  in  column  1).  Table  1 out- 
lines the  content  of  the  RCINPT  data  records.  Figure  7 presents  a schematic 
of  the  deck  setup. 

The  number  of  cards  in  records  1 and  2 and  the  number  of  strings  in  rec- 
ords after  record  2 are  limited  by  core  storage  allotments.  Program  termina- 
tion or  improper  functioning  will  result  if  the  limits  are  exceeded.  The  first 
record  of  data  input  to  the  program  is  a title  card  followed  by  a list  of 
street  names  and  identification  numbers.  The  cards  must  be  placed  in  numeri- 
cal order,  but  need  not  be  in  alphabetical  order,  and  the  numbering  need  not 
be  continuous. 

The  second  data  record  contains  information  that  causes  the  program  to 
generate  from  one  to  ten  maps  that  can  be  used  to  debug  the  map-description 
strings.  Errors  in  the  strings  are  found  by  comparing  the  computer  map  with 
the  original  maps.  The  computer  maps  can  be  drawn  to  a specific  size  so  that 
they  may  be  overlaid  on  the  original  maps.  One  data  card  is  used  to  specify 
each  computer-drawn  map. 

The  third  data  record  describes  the  overall  map  and  its  coordinate  sys- 
tem. Additional  maps  are  described  in  additional  records,  one  record  per  map. 


Reference 

1.  Iuzzolino,  H.  J.,  Air  Foroe  Refuse-Collection  Scheduling  Program, 

CEEDO-TR-77-32 , Tyndall  Air  Force  Base,  Florida,  January  1978. 
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TABLE  1.  RCINPT  DATA  CARDS 


Record 

Card 

. . 

Columns 

. 

Format 

Contents 

1 

i 

1-80 

8A10 

Title 

2 

1-5 

15 

Street  number  (right  justified) 

11-60 

5A1 0 

Street  name  (left  justified) 

The 

above 

card  may 

be  repe 

ated  up  to  300  times. 

End  of 

Recorc 

7-8-9  multipunched  in  column  1. 

2 

1 

1-10 

F10.0 

Minimum  x-coordinate  of  map 

11-20 

F10.0 

Maximum  x-coordinate  of  map 

1 

21-30 

F10.0 

Length  of  map  in  x-direction,  in  inches 

1 

31-40 

F10.0 

Minimum  y-coordinate  of  map 

41-50 

F10.0 

Maximum  y-coordinate  of  map 

51-60 

FI  0.0 

Length  of  map  in  y-direction,  in  inches 

61-70 

F10.0 

Height  of  drum  plot  strips,  in  inches 

71-72 

12 

Sequence  number  of  the  coordinate  system 

The 

above 

card  may 

be  repe 

ated  up  to  10  times. 

End  of 

Record 

7-8-9  multipunched  in  column  1. 

3 

1 

1-10 

FI  0. 6 

Length,  in  inches,  of  street  length  measure- 

ment  unit 

11-20 

F10. 6 

Map  scale,  in  feet/inches 

21-30 

F10.6 

Length  of  map  coordinate  unit,  in  inches 

31-40 

A10 

Coordinate-use  mode  (left  justified) 

2 

1-5 

F5.0 

Default  speed  limit,  in  mph 

6-10 

15 

Default  number  of  directions  streets  can  be 

traveled 

11-15 

15 

Default  number  of  sides  collected  on  one  pass 

16-20 

F5.0 

Default  refuse  quantity  adjustment  factor 

3 

1-80 

Free 

Map-description  strings 

The 

above 

card  may 

be  repe 

;ated  until  500  street  segments  have 

been  described. 

i l 

End  of 

> 

Record 

1 

7-8-9  multipunched  in  column  1. 

i i 

The  above  record  m 
i i 

ay  be  re 

ipeated  for  each  additional  map. 

54 


Figure  7.  RCINPT  Data  Deck  Setup 
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The  first  card  in  each  record  describing  a map  contains  the  unit  (in  inches) 
for  measuring  street  lengths,  the  map  scale  in  feet  per  inch,  the  length  of 
the  MCI)  in  inches,  and  the  mode  of  coordinate  use.  If  street  lengths  can  be 
measured  in  units  precisely  equivalent  to  0.01  mile,  the  unit  of  length  must 
be  left  blank.  The  mode  of  coordinate  use  is  either  AVERAGE  (or  AVG),  FIRST, 
or  LAST.  LAST  should  oe  used  for  debugging  the  maps,  and  AVERAGE  or  FIRST 
should  be  used  for  final  scheduling. 

The  second  data  card  in  each  record  of  map  data  allows  four  default  values 
to  be  redefined.  The  default  values  are  15-mph  speed  limit,  two-way  street, 
two  sides  collected  on  one  pass,  and  an  RQAF  of  1.0.  These  four  pieces  of  in- 
formation-speed limit,  number  of  ways  of  travel,  number  of  sides  collected  on 
one  pass,  and  RQAF--are  needed  for  each  street;  if  they  are  not  specified  on 
the  street-description  cards,  the  default  values  will  be  used.  If  the  defaults 
are  to  be  used,  the  second  card  may  be  blank,  but  it  must  not  be  omitted. 

The  remaining  cards  in  record  3 describe  the  map.  The  description  is  given 
in  pieces  called  strings;  each  string  may  describe  all  or  part  of  a street  and 
thus  may  include  one  or  more  segments  (the  portion  of  a street  between  two 
nodes).  A string  may  occupy  part  of  one  card  or  extend  over  several  cards;  a 
plus  sign  in  column  80  indicates  that  the  string  extends  to  the  next  card.  The 
cards  are  read  in  free  format;  i.e.,  the  numbers  need  not  be  in  specific 
columns.  The  order  is  important  and,  under  certain  circumstances , missing 
items  must  be  indicated  by  commas. 


The  following  items  make  up  a string; 


street  identification  number 

starting  node  number 

length  to  next  node 

refuse  quantity  from  right  side 

refuse  quantity  from  left  side 

next  node  number 

/ 

speed  limit 

number  of  ways  of  travel 

number  of  sides  collected  on  one  pass 

refuse  quantity  adjustment  factor 

( 

x-coordinate  of  starting  node 


This  group  may  be  repeated 
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y-coordinate  of  starting  node 

) 

shape  code 

( 

x-coordinate  of  terminal  node 
» 

y-coordinate  of  terminal  node 

) 

The  first  item  in  a string  is  the  street  identification  number.  The  node  num- 
ber of  an  intersection  or  street  end-point  follows.  The  next  four  items  are 
the  length  to  the  next  node  (in  street-length  units);  the  refuse  quantities  on 
the  right  and  left  sides  of  the  street,  respectively;  and  the  next  node  number. 
Terminal  decimal  points  for  lengths  are  optional.  The  refuse  quantity  and  the 
node  numbers  must  be  positive  integers. 

The  four  items  discussed  in  the  preceding  paragraph  are  the  properties  of 
a segment  and  may  be  repeated  to  add  additional  segments  to  the  string  as  long 
as  the  street  number  initiating  the  string,  the  shape  code,  and  the  four  items 
following  the  slash  (see  next  paragraph)  continue  to  apply  to  the  string  as  a 
whole.  If  a string  must  be  ended  before  all  of  a street  with  a given  name  is 
finished,  the  description  of  that  street  is  resumed  on  the  next  string.  A 
slash  is  punched  after  the  last  node  number  on  a string  unless  all  of  the  re- 
maining items  are  to  be  omitted. 

Following  the  slash  are  the  speed  limit,  number  of  ways  of  travel,  number 
of  sides  serviced  on  one  pass,  and  RQAF.  If  any  of  these  items  are  omitted, 
the  default  values  discussed  earlier  are  used.  Because  these  four  are  posi- 
tional parameters,  missing  preceding  parameters  must  be  indicated  by  commas. 
(Missing  trailing  parameters  need  not  be  indicated  by  commas.)  Since  the  posi- 
tional parameters  must  apply  to  all  segments  of  a given  string,  the  string  must 
end  if  the  next  segment  will  require  a different  value  for  one  of  the  four  pa- 
rameters. If  a street  is  one-way,  the  numeral  1 follows  the  speed  limit  (or 
a comma  indicating  default  speed  limit),  and  the  order  of  nodes  on  the  string 
must  follow  the  direction  of  travel.  Very  wide  streets  and  those  having 
medians  usually  require  refuse  collection  on  only  the  right  side  of  the  street, 
so  a 1 would  be  entered  as  the  third  parameter  after  the  slash.  The  RQAF  is 
used  to  adjust  (by  multiplication)  the  refuse  quantities  for  all  segments  on 
the  string.  For  example,  if  refuse  quantity  has  been  set  at  one  unit  per 
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family,  then  an  RQAF  of  0.5  could  be  used  as  an  approximate  adjustment  for 
areas  housing  one  person  per  dwelling. 

The  coordinates  of  the  first  node  of  the  string,  measured  in  the  coor- 
dinate system  appended  to  the  map,  follow  the  last  parameter  specified  after 
the  slash.  The  two  coordinates  are  separated  by  a comma,  and  the  pair  is  en- 
closed in  parentheses.  No  node  may  have  coordinates  (0.,  0.)  because  (0.,  0.) 
or  (,)  are  used  to  indicate  that  the  coordinates  have  been  specified  previously. 
If  all  four  values  after  the  slash  are  default  values,  the  coordinates  follow 
immediately  after  the  slash. 

If  the  string  does  not  describe  a straight  street,  a shape  code  follows 
the  first  node  coordinates  [or  follows  (,)  if  the  first  node  coordinates  have 
been  given  by  a previous  string,  either  by  explicit  specification  or  by  inter- 
polation]. The  coordinates  of  the  terminal  lode  of  the  string  follow  the 
shape  code.  If  the  coordinates  have  already  been  specified,  nothing  need 
follow  the  shape  code. 

Sample  RCINPT  card  data  for  Kirtland  Air  Force  Base  (Fast)  are  shown  in 
Appendix  D.  Since  only  one  map  was  used,  the  data  consist  of  only  three  rec- 
ords. Three  output  maps  are  produced  by  the  data  cards  in  record  2:  one 
overall  map  and  additional  maps  of  two  housing  areas. 


2.  OUTPUT 

Output  from  program  RCINPT  consists  of  three  disk  files,  one  plot  file, 
and  printed  output.  The  disk  files  are  cataloged  for  later  use  with  programs 
PHASE2,  PHASE3,  and  PHASE4.  The  plot  file  is  either  a disk  file  or  a physical 
tape  depending  on  the  local  system.  The  plot  file  produces  maps  used  to  debug 
input  data.  The  printed  output  is  used  to  debug  the  input  data  and  for  de- 
scriptions of  nodes  and  segments  referred  to  in  later  programs. 


58 


a. 


Disk  Files 


Disk  file  TAPE1  contains  segment  data  and  the  map  distance  conversion 
factor  (miles  per  MCU)  for  the  overall  map.  All  of  the  data  are  written  by  one 
binary  WRITE  statement.  The  first  word  is  the  count  of  the  number  of  segments. 
The  segment  data  follow,  11  words  per  segment  for  each  segment.  After  the  seg- 
ment data  comes  the  overall  distance  conversion  factor.  The  list  used  in  the 
WRITE  statement  is  as  follows: 

KF,  ( ( I ST G ( I , J ) , 1=1 ,11),  J=1 ,KF) , SVAV(l). 

The  ISTG  array  is  equivalenced  to  the  STG  array,  so  some  of  the  contents  in  the 
ISTG  list  are  floating-point.  The  items  in  the  list  are  described  as  follows: 

KF  = number  of  segments 

I STG ( 1 , J ) = street  number 

ISTG(2,J)  = starting  node  number 
ISTG(3,J)  = ending  node  number 
STG(4,J)  = street  length,  in  miles 
ISTG(5,J)  = number  of  houses 
STG(6,J)  = speed  limit,  in  mph 
ISTG(7,J)  = number  of  ways  of  travel 
STG(8,J)  = refuse  quantity  adjustment  factor 
STG(9,J)  = x-coordinate  of  segment  midpoint 
STG(10,J)  = y-coordinate  of  segment  midpoint 

I STG (11  ,J)  = shape  code 

SVAV(l)  = map  distance  conversion  factor,  in 

The  index  J corresponds  to  the  Jth  segment  and  is  also  the  segment  number, 
since  the  program  numbers  the  segments  sequentially,  starting  at  1.  The  sign 
of  the  number  of  refuse  units  is  used  to  indicate  whether  collection  is  on 
both  sides  or  on  only  the  right  side  of  the  segment.  ISTG(5,J)  is  negative 
to  indicate  collection  on  only  the  right  side  of  the  street,  or  positive  to 
indicate  collection  from  both  sides  of  the  street  corresponding  to  segment 
J. 

Disk  file  TAPE2  contains  refuse  quantity  information  and  the  node 
data.  All  of  the  data  are  written  by  one  binary  WRITE  statement.  The  first 
three  words  are  the  total  number  of  houses  or  stops,  the  total  refuse  quantity, 
and  a count  of  the  number  of  nodes.  The  node  data  follow,  four  words  per  node 
for  each  node.  The  list  used  in  the  WRITE  statement  is: 

NHTOT,  TOTREF,  KNODES,  (N0DNUM( I ) ,NBS( I ) ,XN0D( I ) , YN0D( I ) , 1=1,  KNODES). 
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The  items  in  the  list  are  described  as  follows: 

total  number  of  houses  or  stops 
total  refuse  quantity 
number  of  nodes 
node  number 

packed  bounding  segment  numbers  _ These  are  repeated 

x-coordinate  of  node  for  each  node, 

y-coordinate  of  node 

The  index  I corresponds  to  the  Ith  node.  Variable  NBS(I)  contains  up  to  6 
segment  numbers,  each  occupying  10  of  the  60  bits,  for  segments  bounding  node 
NODNUM(I).  Since  node  numbers  are  assigned  by  the  user,  NODNUM(I)  usually  is 
not  the  same  as  I. 

Disk  file  TAPE3  contains  the  street  numbers  and  names  (70  characters 
each).  The  data  are  written  100  streets  at  a time,  using  a parity  1 BUFFER  OUT 
statement.  The  user  should  not  specify  more  than  300  streets.  The  data  writ- 
ten consist  of  all  of  array  NUMSTR  (street  numbers),  followed  by  all  of  array 
NAMSTR  (street  names).  The  output  pointers  are  NUMSTR,  NAMSTR(7,  100).  Since 
the  arrays  are  adjacent  in  storage,  each  record  of  output  appears  as  follows: 

NUMSTR(l)  = number  of  first  street 


NHT0T 

T0TREF 

KN0DES 

N00NUM(I)  = 

NBS(I) 

XNOD(I) 

YNOD(I) 


NUMSTR(IOO)  = number  of  100th  street 
NAMSTR (1 ,1  )~ 

- 7-word  name  of  first  street 

NAMSTR (7,1 ) _ 

NAMSTR(2,1)  = first  word  of  second  street  name 


NAMSTR(7 ,1 00)  = last  word  of  100th  street  name 
The  last  record  contains  zeros  in  unused  words. 

b.  Plot  File 

File  TAPE8,  the  plot  file,  will  be  on  disk  or  tape  depending  on  the 
procedure  used  by  the  local  installation  to  produce  plots.  Each  output  map 
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requested  occupies  one  file.  RCINPT  generates  an  empty  file  to  terminate 
TAPE8.  The  structure  of  the  file  depends  on  the  local  installation. 

c.  Printed  Output 

The  printed  output  consists  of  five  sections:  a listing  of  street 
numbers  and  names,  the  map-description  data  and  associated  error  messages,  a 
summary  of  segment  data,  a summary  of  node  data,  and  some  parameters  used  in 
each  input  and  output  map.  (Appendix  E presents  a sample  of  the  printed  out- 
put.) The  parameters  for  the  input  maps  are  printed  before  the  corresponding 
map-description  data.  The  parameters  for  the  first  output  map  follow  the 
first  map-description  string.  The  parameters  for  the  other  output  maps  are 
printed  after  the  node  data. 

The  street  listing  gives  the  street  numbers  and  names,  one  pair  per 

line. 


In  the  map-description  data  section,  the  parameters  for  input  maps 
are  preceded  by  the  heading  "PARAMETERS  FOR  MAP  n"  (n  is  the  sequence  number 
of  the  map).  Eight  of  the  next  ten  lines  give  the  values  determined  by  the 
first  two  data  cards  of  the  map-description  data  records.  Two  additional 
pieces  of  information  are  given:  the  coordinate  scale  factor  in  miles  per 
MCI)  and  the  ratio  of  the  current  scale  factor  to  that  for  the  overall  map. 

When  fields  of  the  second  data  card  of  record  3 are  blank,  program-defined 
default  values  are  used  and  are  printed  with  the  map  parameters. 

Seventeen  map  parameters  are  printed  for  output  maps.  SCR  is  the 
scale  ratio.  It  should  be  1.0  for  the  first  output  map  and  usually  will  be 
between  0.04  and  1.0  for  subsequent  maps.  TX  and  TY  are  the  x-  and  y-components 
of  translation  of  the  current  coordinate  system  with  respect  to  the  overall 
(first)  coordinate  system.  These  are  measured  in  MCU  in  the  overall  system. 

The  minimum  and  maximum  x-  and  y-coordinates  are  given  for  the  output  map, 
first  as  XMIN,  XMAX,  YMIN,  and  YMAX,  measured  in  MCU  for  the  coordinate  system 
specified  in  column  72  of  the  corresponding  mao  output  card,  and  again  as  XL, 
XR,  YB,  and  YT,  measured  in  MCU  in  the  overall  coordinate  system.  Variables 
XSC  and  YSC  identify  the  x and  y map  scales  in  inches  per  MCU.  The  plot  strip 


height  (PHGT)  and  total  length  (PLEN)  are  printed  in  inches.  The  height  of 
the  plot  or  the  height  of  a strip,  whichever  is  smaller,  is  printed  as  YCUT, 
in  MCU.  The  last  item  is  CNVLEN,  the  length  conversion  in  miles  per  street- 
length  unit.  Usually  CNVLEN  is  between  0.001  and  0.02.  As  each  output  map 
after  the  first  is  completed,  the  message  "COMPLETED  PLOT  OF  MAP  n"  is  printed 
(n  is  the  sequence  number  of  the  map).  No  message  is  printed  when  the  first 
map  is  completed. 

The  map-description  cards  are  printed  exactly  as  they  are  read. 
Sample  output  is  shown  on  page  2 of  Appendix  E.  Error  messages  are  printed 
immediately  following  any  map-description  string  that  causes  any  problem  for 
program  RCINPT.  These  messages  are  described  in  Section  VII. 

The  segment  data  summary  printout  starts  with  the  problem  title, 
the  number  of  houses  serviced,  the  total  refuse  quantity,  and  a count  of  the 
segments.  The  word  HOUSES  may  sometimes  actually  indicate  the  number  of 
stops,  the  number  of  families  serviced,  the  number  of  containers  serviced, 
the  refuse  volume,  or  the  refuse  weight,  depending  on  the  meaning  of  the 
refuse  data  on  the  map-description  strings.  The  total  refuse  quantity  is  the 
sum  of  RQAF  times  the  refuse  numbers  on  the  map-description  strings.  The  seg- 
ment data  consist  of  the  11  items  per  segment  written  to  disk  file  1,  preceded 
by  a segment  sequence  number  that  hereafter  will  identify  that  particular  seg- 
ment. The  shape  code,  in  the  column  headed  SF,  is  printed  in  octal  because  it 
may  be  either  two  characters  (16  zeros  followed  by  4 octal  digits)  or  a 
floating-point  number  (20  octal  digits  beginning  with  a 1 or  a 6).  A state- 
ment at  the  end  of  the  segment  list  indicates  that  if  a negative  number  of 
houses  is  entered  in  the  column  headed  NH,  collection  is  from  the  right  side 
of  the  street,  only. 

The  node  list  begins  with  a problem  title.  The  columns  are  headed 
I, NODE, T,  T*X,  T*Y,  NBR1 , NBR2,  NBR3,  NBR4,  NBR5,  and  NBR6.  The  column  headed 
I is  the  line  number,  and  the  column  headed  NODE  is  the  node  number.  The 
column  headed  T gives  the  number  of  times  the  node  was  referenced.  X*T  is  the 
sum  of  the  T values  of  the  node's  x-coordinate.  Y*T  is  the  corresponding  value 
for  the  y-coordinate.  NBR1  through  NBR6  are  segment  numbers  for  segments 
bounding  the  node.  Zeros  are  printed  when  there  are  fewer  than  six  neighbor- 
ing segments. 
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SECTION  V 

PROGRAM  REQUIREMENTS 


1.  SYSTEM 

Program  RCINPT  is  written  entirely  in  FORTRAN  IV.  The  program  runs  on  a 
CDC  6600  using  a SCOPE  3.4.4  operating  system. 

Eleven  obvious  types  of  computer-dependent  coding  are  used  in  program 
RCINPT  and  its  subroutines.  A 132-character  output  line  is  assumed  in  the 
main  program.  Ten  characters  per  word  are  assumed  in  the  main  program  and  in 
subroutines  STRINP  and  NUMBER.  A 60-bit  word  is  assumed  in  the  main  program. 
System  function  SHIFT  is  used  by  the  main  program  and  by  subroutines  SHAPCOM, 
LINEIN,  AXIS,  and  NUMBER.  An  R-format  is  used  in  subroutines  SHAPCOM,  LINE  IN , 
and  COORD.  A non-zero  floating-point  characteristic  is  assumed  in  subroutine 
SHAPCOM.  Six-bit  characters  are  assumed  in  subroutines  SHAPCOM,  LINEIN,  AXIS, 
and  NUMBER.  Further,  the  CDC  values  for  the  characters  are  assumed  in  subrou- 
tine AXIS.  Multiple  entry  points  are  used  in  subroutine  LINEIN.  An  ENCODE 
statement  is  used  in  subroutine  NUMBER.  A dollar  sign  is  used  between  FORTRAN 
statements  in  program  RCINPT  and  subroutines  STRINP,  MAPGRID,  AXIS,  MAPPLT, 
SHAPCOM,  and  COORD. 

More  subtle  types  of  machine  dependencies  may  exist,  according  to  the 
machine  used. 

2.  STORAGE 

The  core  requirement  is  slightly  less  than  72,0008  words,  but  may  vary 
somewhat  in  accordance  with  the  plotting  system  used  by  the  local  installa- 
tion. The  peripheral  storage  for  output  disk  files  should  not  exceed  9,000 
words  for  TAPE1,  2,600  words  for  file  TAPE2,  1,000  words  for  file  TAPE3,  and 
1.5  million  words  for  file  TAPE8,  the  plot  file.  The  plot  file  will  usually 
contain  about  20,000  words  per  output  map. 
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3. 


TIME 


At  least  70  percent  of  the  program  execution  time  is  used  to  generate 
the  plot  file.  The  execution  CP  time  is  approximately  .02  STnT  + .04  2 $. 

seconds,  where  S^-gy  is  the  total  number  of  segments  and  S.  is  the  number  of 

segments  on  the  i^h  output  map.  It  is  time-consuming  to  count  the  total  num- 
ber of  segments  and  the  number  of  segments  on  each  map;  therefore,  60  seconds 

can  be  used  as  a reasonable  upper  limit.  The  maximum  possible  CP  time  using 
10  full  maps  of  a 720-segment  network  would  be  302.4  seconds,  but  it  is  not 
likely  that  the  maximum  will  be  needed  because  9 of  the  10  maps  usually  will 
be  of  nonoverlapping  regions.  Upper  limits  on  10  and  PP  time  are  roughly 
60  and  100  seconds,  respectively.  Not  enough  data  are  available  to  determine 
even  a rough  functional  relation  for  10  and  PP  time. 


64 


SECTION  VI 
PROGRAM  LIMITATIONS 


1.  PLOTTER 

The  plotter  used  by  RCINPT  can  be  a drum  plotter  with  at  least  a 10-inch 
drum  or  a flatbed  plotter  with  a 30-inch-square  or  larger  bed.  The  program 
assumes  that  a 30-inch  height  is  available  if  the  user  does  not  indicate  other- 
wise on  the  output  description  cards.  Maps  may  not  be  drawn  in  strips  on  a 
flatbed  plotter;  therefore,  the  height  and  width  specified  for  the  output  map 
must  not  be  greater  than  the  plot  strip  height  (the  7th  field  on  output  map- 
description  cards).  Some  systems  limit  plot  lengths;  for  example,  Calcomp 
plots  may  not  be  longer  than  120  inches  at  the  Air  Force  Weapons  Laboratory. 
Each  strip  of  an  output  map  has  a 1-inch  space  after  it,  with  an  additional 
inch  after  the  last  strip.  Thus  a 30-  by  30-inch  map  plotted  as  three  30- 
by  10-inch  strips  generates  a plot  94  inches  long.  The  user  must  consider 
any  limitations  on  plot  lengths  imposed  by  the  local  system  when  the  output 
map-description  cards  are  punched. 


2.  SEGMENTS 

Program  RCINPT  can  read  map-description  strings  containing  up  to  720  seg- 
ments. Beyond  this  point,  the  program  will  attempt  to  store  data  beyond  its 
field  length,  causing  program  termination  and  a M0DE1  error  (address  out  of 
range).  Each  segment  that  requires  collection  from  separate  sides  on  separate 
passes  counts  as  two  segments.  If  more  than  500  segments  are  used,  the  map 
description  will  have  to  be  simplified  because  program  PHASE3  will  accept  only 
500  segments.  If  segment  storage  is  increased  in  programs  PHASE2,  PHASE3,  and 
PHASE4,  a limit  of  1023  segments  is  imposed  by  the  packing  used  in  variable 
NBS  (each  segment  number  is  allowed  to  fill  10  bits).  No  more  than  six  seg- 
ments may  meet  at  one  node. 
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3. 


NODES 


Program  RCINPT  has  storage  set  aside  for  500  nodes.  If  more  than  500 
nodes  are  read,  output  map-description  data  will  be  overwritten.  This  situa- 
tion usually  causes  the  program  to  terminate  while  the  first  output  map  is 
being  plotted.  If  the  program  completes  the  first  map,  the  node  count  at  the 
beginning  at  the  node  listing  will  indicate  the  excessive  number  of  nodes. 

Any  maps  drawn  will  be  unreliable.  No  other  message  is  given.  The  map  de- 
scription will  have  to  be  simplified,  or  else  node  storage  must  be  increased 
in  the  arrays  in  COMMON  block  NDDATA. 

Node  numbers  must  be  between  1 and  99999,  inclusive.  Numbering  need  not 
be  consecutive.  The  five-digit  limitation  comes  from  the  format  used  to  print 
the  node  listing.  In  addition  to  this  limitation,  subroutine  LINEIN  limits 
the  node  number  to  15  digits  by  the  algorithm  used  to  build  integers.  Non- 
positive node  numbers  will  cause  string  scanning  to  terminate,  and  various 
unpredictable  errors  will  result. 


4.  STREETS 

Street  number  and  name  data  should  be  limited  to  300  cards.  The  numbers 
may  be  between  1 and  99999,  inclusive;  the  names  may  be  70  characters  long, 
although  only  the  first  30  characters  are  used  by  PHASE4  on  maps  and  printed 
schedules.  The  cards  must  be  in  increasing  numerical  order,  but  not  neces- 
sarily in  alphabetical  order.  The  numbers  need  not  be  consecuti ve--gaps  are 
al lowed. 

If  more  than  300  streets  are  used,  the  excess  will  be  ignored  by  program 
PHASE4.  No  warning  is  given  by  program  RCINPT.  If  street  numbers  are  not  in 
order,  some  street  names  will  be  missing  from  the  printed  output  and  map  pro- 
duced by  PHASE4,  but  no  warning  is  given. 

A street  name  may  have  more  than  one  number.  This  will  cause  the  name  to 
be  appended  to  PHASE4  route  maps  once  for  each  number.  Each  street  number  must 
occur  only  once;  otherwise  all  but  one  of  the  occurrences,  not  necessarily  the 
first,  will  be  ignored  by  PHASE4. 
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5.  OUTPUT  MAPS 


1 


. The  number  of  output  maps  must  be  between  1 and  10,  inclusive.  If  no 

maps  are  specified,  the  program  will  terminate  in  subroutine  MAPPLT  while  try- 
ing to  use  the  result  of  a division  by  zero.  If  more  than  10  maps  are  speci- 
fied, the  program  will  terminate  either  from  an  end-of-record  encountered  af- 
ter the  11th  card  or  from  an  improper  format  on  the  12th  card.  No  specific 
message  is  printed  if  the  end-of-record  is  encountered,  but  the  program  will 
print  segment  and  node  counts  of  zero. 


6.  ANGLE  SHAPE  CODE 

Angle  shape  codes  are  limited  to  eight  characters.  Any  longer  shape 
code  will  be  misinterpreted  by  subroutine  SHAPCOM.  No  warning  is  printed. 
This  restriction  should  not  be  bypassed  because  the  first  12  bits  of  each 
angle  shape  code  are  tested  by  SHAPCOM  to  determine  whether  the  shape  code 
is  in  character  or  floating-point  form.  Also,  angle  shape  codes  of  five  or 
six  characters  (letter  and  decimal  point  included)  would  be  more  precise 
than  any  other  measurement,  so  more  characters  are  not  necessary. 


| (The  reverse  of  this  page  is  blank.) 
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SECTION  VII 

ERROR  MESSAGES  AND  CORRECTIVE  ACTION 


1.  A PARITY  ERROR  OCCURRED  DURING  THE  BUFFER  OUT  TO  UNIT  3 OF  THE 
PRECEDING  nnn  STREETS. 


Note: 

nnn  will  be  some  number  from  1 to  100. 

Type: 

Warning. 

Source: 

Subroutine  STRINP. 

Location: 

In  street  listing. 

Meaning: 

Self-explanatory. 

Action: 

Either  rerun  the  job  or  use  the  street-name  cards,  rather  than 
TAPE3,  as  input  to  program  PHASE4. 

2.  NO  STREET  NUMBERS  AND  NAMES  WERE  SPECIFIED. 


Type: 
Source: 
Location: 
Meani ng: 
Action: 


Warning. 

Subroutine  STRINP. 

Where  street  listing  would  normally  appear. 

Self-explanatory. 

If  street  names  are  intentionally  absent,  map-description  strings 
must  not  include  street  numbers.  If  street  names  are  included  in 
the  data  deck,  check  for  improper  location  of  these  cards  or  for 
an  extra  end-of-record  card. 


3.  ---ONLY  ONE  MAP  IS  ALLOWED  WHEN  THE  VARIABLE  MAP  COORDINATE 
OPTION  IS  USED. 

JOB  TERMINATED. 

Type:  Fatal. 

Source:  Program  RCINPT. 

Location:  In  the  map-description  string  description. 

Meaning:  If  the  map  scale  is  blank,  one  map,  portions  of  which  are  drawn  to 

different  scales,  may  be  used  (this  option  is  untested  and  not 
recommended).  Only  one  record  of  map  descriptions  may  then  be  used; 
the  message  indicates  that  a second  record  has  been  encountered. 

In  other  cases,  the  map  scale  has  been  omitted  or  mispunched. 
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Action: 


If  the  variable  map-scale  option  is  desired,  use  only  one  map- 
description  record  followed  by  two  end-of-record  cards  or  an  end- 
of-file  card.  Otherwise,  determine  whether  the  map  scale  has  been 
punched  properly. 


4.  PROBLEMS  IN  ABOVE  CARD. 


Type: 

Source: 

Location: 

Meaning: 


Action: 


Fatal;  job  terminates  after  20  such  errors  or  when  all  map- 
description  strings  have  been  processed,  whichever  occurs  sooner. 
Program  RCINPT. 

Following  the  printing  of  the  map-description  string. 

The  form  of  the  preceding  map-description  string  is  improper.  Pos- 
sibilities include  missing  or  extra  fields,  a zero  node  number,  mis- 
punched  fields,  or  illegal  characters. 

Check  the  string  against  the  form  given  in  Section  IV,  item  by  item. 
Also  check  for  a plus  sign  in  column  80  if  the  string  is  to  be  con- 
tinued, or  for  an  end-of-record  where  a continuation  card  is  ex- 
pected. If  these  procedures  do  not  disclose  the  error,  find  vari- 
able II  in  a dump  of  LINEIN  to  locate  the  column  before  which  the 
problem  occurred.  Also  find  variable  IBRK  in  RCINPT  to  get  an  idea 
of  the  error  type.  If  IBRK  > 0,  a zero  node  number  is  indicated. 

IBRK  = 0,  which  should  be  impossible,  indicates  that  LINEIN  resumed 
the  scan  of  a card  after  column  80.  IBRK  = -1  shows  that  an  end-of- 
record  has  been  encountered  unexpectedly.  IBRK  = -2  indicates  that 
a number  was  not  in  the  form  expected  by  LINEIN. 


5.  BREAK  = a IN  THE  ABOVE  LINE  FOUND  WHILE  SCANNING  FOR  ( BEFORE 
FIRST  COORDINATE. 


Note: 

Type: 

Source: 

Location: 

Meaning: 


a may  be  any  CDC  character  other  than  (. 

Warning,  but  the  remainder  of  the  card  is  ignored. 

Program  RCINPT. 

Following  the  printing  of  the  map-description  string. 

The  first  node's  coordinates  may  be  missing,  or  an  extra  character 
may  have  occurred  before  the  first  node's  coordinates. 
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Action:  Append  the  first  node's  coordinates  or  (,)  if  the  node  occurred  on  a 

previous  string.  If  the  node's  coordinates  are  already  present,  the 
order  of  the  items  on  the  string  may  be  incorrect  or  an  item  may  have 
been  mispunched. 


6.  ---THE  ABOVE  STRING  DOES  NOT  START  OR  END  ON  A PREVIOUSLY 
DEFINED  NODE. 


Type: 

Source: 

Location: 

Meaning: 


Action: 


Fatal;  job  terminates  after  20  such  errors  or  when  all  map- 
description  strings  are  processed,  whichever  occurs  sooner. 

Program  RCINPT. 

Following  the  printing  of  the  first  map-description  string  of  any 
input  map  other  than  the  first. 

The  first  string  of  the  record  starts  and  ends  on  nodes  not  yet  de- 
fined to  the  program.  One  node  is  needed  to  obtain  the  translation 
between  the  current  map-coordinate  system  and  the  first  map- 
coordinate  system. 

If  neither  the  starting  nor  the  ending  node  has  been  used  previously, 
either  add  one  of  them  to  a previous  map  description  or  start  the 
current  map  description  with  a string  that  starts  or  ends  on  a node 
used  previously.  If  the  starting  or  ending  node  has  been  used  pre- 
viously, an  error  in  that  string  has  prevented  the  correct  processing 
of  the  coordinates,  and  that  string  must  be  corrected.  Coordinates 
of  (0.,  0.)  must  not  be  used  for  any  node. 


7.  NO  COORDINATES  WERE  GIVEN  FOR  NODE  nnnnn.  THEY  WILL  BE  ASSUMED 
TO  BE  (0,0). 


Note: 

Type: 

Source: 

Location: 

Meaning: 

Action: 


nnnnn  will  be  a number  from  1 to  99999. 

Warning,  but  the  output  map  will  be  incorrect. 

Program  RCINPT. 

Following  the  printing  of  a map-description  string. 
Self-explanatory. 

Coordinates  must  be  specified  for  the  node  unless  an  error  in  a 
previous  string  containing  the  node  caused  the  problem. 


8.  ---BAD  DISTANCE  SPECIFICATION  ON  PREVIOUS  LINE. 

THE  MAP  COORDINATE  SCALE  (dd.ddd)  DEVIATES  FROM  THE  DEFAULT 
VALUE  (dd.ddd  MILES  PER  MAP  COORDINATE  UNIT). 

I 

Note:  dd.ddd  represents  a floating-point  number. 

Type:  Warning. 

Source:  Program  RCINPT. 

Location:  Following  the  printing  of  a map-description  string. 

Meaning:  The  distance  scale  implied  by  length  in  the  string  differs  appre- 

ciably from  the  scale  in  use  for  the  rest  of  the  map  because  either 
length  or  coordinates  were  incorrectly  specified.  (In  some  cases, 
the  incorrect  specification  may  have  been  given  deliberately;  it 
may  be  desirable  to  have  a given  street  drawn  other  than  to  scale.) 

Action:  Errors  in  length  that  cause  scale  changes  of  less  than  20  percent 

can  usually  be  tolerated.  Small  errors  in  coordinate  specifications 
can  be  tolerated.  If  the  output  map  looks  excessively  distorted  or 
if  very  precise  distance  measurements  are  needed,  the  lengths  and 
coordinates  in  the  string  should  be  corrected. 

9.  THE  PREVIOUS  LINE  CHANGES  THE  MAP  COORDINATE  SCALE  FROM  dd.ddd 
TO  dd.ddd  MILES  PER  MAP  COORDINATE  UNIT. 

Note:  dd.ddd  represents  a floating-point  number. 

Type:  Warning. 

Source:  Program  RCINPT. 

Location:  Following  the  printing  of  a map-description  string. 

Meaning:  The  straight  string  indicated  has  a total  length  or  coordinates 

that  cause  a significant  change  in  the  map  scale.  The  new  scale 
will  be  used  until  another  significant  change  is  caused  by  another 
straight  string. 

Action:  If  the  scale  change  is  intentional,  no  action  is  required;  if 

unintentional,  the  lengths  or  coordinates  on  the  string  must  be 
corrected. 

10.  ---BAD  DISTANCE  SPECIFICATION  IN  PREVIOUS  LINE . — 

MAP  DISTANCE  (dd.ddd)  EXCEEDS  TOTAL  SEGMENT  LENGTH  (dd.ddd) 
(nnn.nn,nnn.nn)  TO  (nnn.nn,nnn.nn) 

Notes:  dd.ddd  represents  floating-point  distances  in  miles;  nnn.nn  repre- 

sents floating-point  node  coordinates  in  MCU. 
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Type:  Warning. 

Source:  Subroutine  SHAPCOM. 

Location:  Following  the  printing  of  a map-description  string  or  following  the 
parameters  for  some  output  map  after  the  first. 

Meaning:  The  string  is  shorter  than  the  distance  between  the  end-point  nodes. 
Action:  Small  differences,  perhaps  0.02  mile  or  less,  usually  are  tolerable. 

Unacceptable  differences  require  a correction  of  length  or  coordin- 
ates. 


11.  ILLEGAL  CHARACTER  a IN  SHAPE  FACTOR  IN  ABOVE  LINE. 

Note:  a is  any  CDC  character. 

Type:  Warning. 

Source:  Subroutine  SHAPCOM. 

Location:  Following  the  printing  of  a map-description  string. 

Meaning:  Self-explanatory. 

Action:  The  program  will  treat  the  string  as  a straight  string.  The  shape 

code  should  be  corrected. 

12.  BAD  ANGLE.  SIDES  = dd.ddd  dd.ddd  SPAN  = dd.ddd 
THE  ANGLE  WILL  BE  TREATED  AS  A STRAIGHT  LINE. 

Note:  dd.ddd  represents  floating-point  lengths  in  miles. 

Type:  Warning. 

Source:  Subroutine  SHAPCOM. 

Location:  Following  the  printing  of  a map-description  string  with  an  angle 
shape  code. 

Meaning:  The  triangle  inequality  has  been  violated;  specifically,  the  total 

length  of  two  sides  of  the  triangle  formed  by  connecting  the  end 
points  of  the  angle  is  less  than  the  length  of  the  third  side,  which 
is  geometrically  impossible. 

Action:  The  lengths  and  coordinates  of  the  string  should  be  checked  and 

corrected. 

13.  ERROR  IN  COLUMN  nn  OF  FOLLOWING  LINE. 

Note:  nn  indicates  a column  number  from  1 to  80.  Column  1 of  the  next  line 

is  preceded  by  an  asterisk. 
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Type: 

Source: 

Location: 

Meaning: 

Action: 


Fatal;  processing  terminates  after  20  such  errors  or  when  map- 
description  processing  is  complete,  whichever  occurs  sooner. 
Subroutine  LINEIN. 

In  the  map-description  string  printout. 

LINEIN  has  encountered  an  illegal  character  in  the  column  indicated 
while  processing  a field. 

The  field  should  be  corrected  if  in  error.  If  not,  a field  may  be 
missing,  or  an  extra  field  may  be  present  earlier  in  the  card. 


14.  STOP  567. 


Type: 

Source: 

Location: 

Meaning: 

Action: 


Fatal . 

Subroutine  AXIS. 

Dayfi le. 

An  illegal  character  has  been  encountered  in  the  format  for  AXIS 
numbers. 

This  field  is  specified  explicitly  in  each  call  to  AXIS  and  is  not 
accessible  to  the  user  through  data  cards.  If  the  field  has  been 
changed  by  the  user,  it  must  start  with  E,  F,  G,  or  I.  A parity 
error  in  the  FORTRAN  source  card  may  have  occurred,  in  which  case 
the  program  should  be  recompiled.  The  format  may  have  been  over- 
written by  some  data  exceeding  its  allocated  storage;  in  this  case, 
the  storage  overflow,  which  will  probably  be  difficult  to  find, 
must  be  corrected. 


Finally,  one  other  error  condition  for  which  no  error  message  is  printed 
can  occur.  If  more  than  500  nodes  have  been  used,  the  output  maps  will  be 
unreliable.  Therefore,  as  a matter  of  routine  the  count  of  nodes  and  segments 
on  the  appropriate  portions  of  the  printed  output  should  be  checked.  If  neces- 
sary, the  map  must  be  simplified  to  place  the  count  within  the  limits  of  500 
segments  and  300  nodes. 
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SECTION  VIII 

RECOMMENDED  PROGRAM  CHANGES 

The  following  changes  in  program  RCINPT  are  recommended: 

1.  Limit  street  names  to  30  characters,  the  maximum  number  PHASE4 
can  use.  This  modification  involves  changing  the  dimension  of 
NAMSTR  to  (3,100)  in  STRINP.  Corresponding  changes  to  coding 
in  STRINP  are  required. 

2.  Add  warning  messages  in  subroutine  STRINP  to  indicate  that 
street  numbers  are  out  of  order  or  that  there  are  more  than 
300  street  name  cards. 

3.  Terminate  reading  of  map-description  strings  in  RCINPT  when 
more  than  700  segments  have  been  read.  Print  a warning  message. 

4.  Terminate  reading  of  map-description  strings  in  RCINPT  when 
more  than  500  nodes  have  been  read.  Print  a warning  message. 

5.  Change  the  variables  that  hold  the  number  of  houses  (NHL, 
ISTG(NH, . . . ) , and  NHT0T)  to  floating-point  so  that  refuse  meas- 
urements with  fractional  parts  need  not  be  scaled  to  integers. 
Similar  changes  will  be  necessary  in  programs  PHASE2,  PHASE3, 
and  PHASE4. 
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END 


Call  PLOT  to  draw  tic  marks 
and  axis. 
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marks. 
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APPENDIX  C 


DEFINITIONS  OF  IMPORTANT  VARIABLES 


Page 

Subrouti ne 

MOVE  5 

136 

Function  IFIND 

136 

Subrouti ne 

STRINP 

136 

Subroutine 

NUMBER 

136 

Subrouti ne 

AXIS 

136 

Subroutine 

MAPGRID 

137 

Subroutine 

SHAPCOM 

137 

Subroutine 

COORD 

138 

Subroutine 

MAPPLT 

139 

Subroutine 

LINEIN 

139 

Program  RCINPT 

140 

135 


Note: 


Al,... 

IF 

II 


IARRAY 

LEN 

NUM 


MSTINC 

NAMSTR 

NS 

NUMSTR 


FORM 

NUM 

TEXT 


ANGAX 

ANGNM 

CAX 

DLBL 


1 

A single  variable  symbol  may  have  different  meanings  in  relation  to  the 
various  subroutines.  For  this  reason,  variables  are  defined  below  for 
each  subroutine  and  for  program  RCINPT. 

SUBROUTINE  M0VE5 

A5  Arrays  to  be  moved 

Subscript  data  are  moved  to 
Subscript  data  come  from 

FUNCTION  IFIND 

Array  being  searched 
Length  of  IARRAY 
Number  being  sought 

SUBROUTINE  STRINP 

Number  of  streets  that  can  be  read  before  writing  data  to  disk 
Array  of  street  names 
Count  of  streets 
Array  of  street  numbers 

SUBROUTINE  NUMBER 

Output  format  for  number 

Number  to  be  plotted 

Character  representation  of  number 

SUBROUTINE  AXIS 

Axis  angle,  in  degrees 
Numbering  angle,  in  degrees 
Cosine  of  axis  angle 
Spacing  of  numbered  tic  marks 


136 


SUBROUTINE  AXIS  (Cont'd.) 


HGT 

Character  height  of  numbering, 

in  inches 

IV 

Value  of  number,  when  integer 

MODE 

Type  of  tic  mark  plotting  indicator 

SAX 

Sine  of  the  axis  angle 

SCALE 

Scale,  in  inches  per  axis  unit 

TIC 

Spacing  of  small  tic  marks 

V 

Value  of  number,  when  floating- 

-point 

VF 

Value  of  end  of  axis 

VI 

Value  of  start  of  axis 

X 

X-coordinate  of  start  of  axis 

Y 

Y-coordinate  of  start  of  axis 

! 

SUBROUTINE  MAPGRID 

NH  Number  of  integral  map  coordinate  units  in  the  vertical  extent 

of  a map  strip 

NPL  Number  of  map  strip 

NV  Number  of  integral  map  coordinate  units  in  the  horizontal  extent 

of  a map  strip 

XDISPL  Displacement,  in  plotter  inches,  from  beginning  of  first  strip 

to  beginning  of  current  map  strip 

XH  X-coordinate,  in  plotter  inches,  of  start  of  a vertical  axis 

YH  Y-coordinate,  in  plotter  inches,  of  start  of  a horizontal  axis 

SUBROUTINE  SHAPCOM 

AVMD  Map  distance  conversion  factor,  in  miles  per  map  coordinate  unit 

BR1  Distance  from  starting  node  to  first  break  in  segment  shape,  in 

miles 

BR2  (Rectangle)  Distance  from  starting  node  to  second  break  in  seg- 

ment shape,  in  miles 

(Angle)  Distance  from  vertex  to  ending  node,  in  miles 

Cll  X-component  of  vector  from  center  to  starting  point  of  S-curve 

or  circular  arc 

C12  Y-component  of  vector  from  center  to  starting  point  of  S-curve 

or  circular  arc 
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SUBROUTINE  SHAPCOM  (Cont'd.) 


CNVLEN 

D 

DD 

H 


ISF 

R 

RPR 

S 

SCR 

SF 

THETA,  6 

TOTLEN 

XNF 

XNI 

YNF 

YNI 


Segment  length  conversion,  in  miles  per  street-length  unit 

Distance,  in  miles,  from  starting  to  stopping  nodes 

Half  the  distance  from  starting  to  stopping  points 

(Circular  Segment)  Distance  from  center  of  circle  to  line  con- 
necting starting  and  stopping  points 

(Angle  Segment)  Distance  from  vertex  to  line  connecting  start- 
ing and  stopping  points 

Shape  code 

Radius  of  curvature  of  circular  segments,  in  miles 
Reciprocal  of  radius  of  curvature 
Perimeter  measurement 

Ratio  of  current  map  distance  conversion  to  that  of  first  map 

Shape  code  for  angles  after  conversions  to  floating  point 

Slope  of  line  from  starting  to  ending  node,  in  radians 

Total  length  of  segment,  in  miles 

X-coordinate  of  ending  node 

X-coordinate  of  starting  node 

Y-coordinate  of  ending  node 

Y-coordinate  of  starting  node 


BR1 

BR2 

CUMLEN 

RPR 

S 

SF 

XNF 

XNI 

YNF 

YNI 


SUBROUTINE  COORD 

Distance  to  first  break  in  segment  shape,  in  miles 

Distance  to  second  break  in  segment  shape,  in  miles 

Cumulative  length  along  string,  in  miles 

Reciprocal  of  radius  of  curvature  of  a circular  segment 

Distance  along  segment 

Shape  code 

X-coordinate  of  ending  node 
X-coordinate  of  starting  node 
Y-coordinate  of  ending  node 
Y-coordinate  of  starting  node 
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SUBROUTINE  MAPPLT 


AVMD 

CNVLEN 

CUMLEN 

INB 

ISF 

ISTG 

KNODES 

MSEQ 

NMAP 

NMAPO 

NODNUM 

NPPSEG 

PHGT 

PLEN 

SCR 

STG 

SVAV 

TOTLEN 

TRX 

TRY 

YCUT 


IAB 

IBRK 

IC 

ICHAR 

ISGN 

ITYP 

IV 

LB 

NT 


Map  distance  conversion,  in  miles  per  map  coordinate  unit 
Street  length  conversion,  in  miles  per  street  length  unit 
Cumulative  street  or  string  length,  in  miles 
Point  within  map-bounds  indicator 
Shape  code 

Array  of  integer-valued  segment  data 
Count  of  nodes 

Sequence  number  of  input  map  coordinate  system 

Map  strip  number  of  current  point 

Map  strip  number  of  previous  point 

Array  of  node  numbers 

Number  of  points  plotted  per  segment 

Height  of  map  strip,  in  inches 

Total  length  of  all  plot  strips,  in  inches 

Ratio  of  current  map  distance  conversion  to  that  of  first  map 

Array  of  floating  point-valued  segment  data 

Array  of  map  distance  conversion  factors 

Total  segment  or  string  length,  in  miles 

Array  of  x-components  of  map  translations  relative  to  overall 
coordinate  system 

Array  of  y-components  of  map  translations  relative  to  overall 
coordinate  system 

Height  of  map  output  strips,  in  map  coordinate  units 


SUBROUTINE  LINEIN 

Break  character  for  bounded  character  string 
Break  character  at  end  of  LINEIN  processing 
Array  of  characters  from  input  line 
Character  currently  being  processed 
Array  of  signs  used  to  build  a number 
Type  of  word  being  processed 
Array  of  values  used  to  build  a number 
Blank  previous  character  indicator 
Number  of  digits  following  a decimal  point 


PROGRAM  RCINPT 


AVMD 

CNVLEN 

CTYPE 

CUMLEN 

FNPT 

INPT 

IQUIT 

ISTG 

KF 

KNODES 

LEN 

MAXERR 

MAXSTG 

MDFILE 

MODE 

NH 

NHTOT 

NN1 

NN2 

NODEF 

NODE  I 

NODNUM 

NRQF 

NSF 

NSPD 

NSTR 

NWAY 

NXMID 

NYMID 

SCALEC 

SCALEM 

STG 


Average  map  distance  conversion  factor,  in  miles  per  map 
coordinate  unit 

Length  conversion  factor,  in  miles  per  street-length  unit 

Mode  of  coordinate  use 

Cumulative  length  along  a string 

Array  for  floating-point  numbers  from  LINEIN 

Array  for  integer  numbers  from  LINEIN 

Count  of  fatal  errors  in  map  description  strings 

Array  of  integer-valued  segment  data 

Count  of  segments 

Count  of  nodes 

Symbolic  subscript  of  STG  array  for  street  length 

Number  of  fatal  errors  allowed  before  program  termination 

Maximum  number  of  segments 

Current  number  of  map-description  record 

String  print  control 

Symbolic  subscript  of  ISTG  array  for  number  of  houses 
Total  number  of  houses 

Symbolic  subscript  of  ISTG  array  for  starting  node  number 

Symbolic  subscript  of  ISTG  array  for  ending  node  number 

End  node  number  of  a string 

Start  node  number  of  a string 
Array  of  node  numbers 

Symbolic  subscript  of  STG  array  for  refuse  quantity  adjustment 
factor 

Symbolic  subscript  of  ISTG  array  for  shape  code 

Symbolic  subscript  of  STG  array  for  speed  limit 

Symbolic  subscript  of  ISTG  array  for  street  number 

Symbolic  subscript  of  ISTG  array  for  number  of  ways  of  travel 

Symbolic  subscript  of  STG  array  for  x-coordinate  of  segment 
midpoint 

Symbolic  subscript  of  STG  array  for  y-coordinate  of  segment 
midpoint 

Coordinate  scale,  in  plotter  inches  per  unit 
Map  scale,  in  true  feet  per  map  inch 
Array  of  floating  point-valued  segment  data 
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PROGRAM  RCINPT  (Cont'd.) 


SVAV  Array  of  map  distance  conversion  factors 

TOTLEN  Total  string  length,  in  miles 

TOTREF  Total  refuse  quantity 

XNOD  Array  of  node  x-coordi nates 

YHCUT  Height  of  output  map  strips,  in  inches 

YNOD  Array  of  node  y-coordi nates 


i 
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APPENDIX  D 
SAMPLE  INPUT  DATA 
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sifter  name 


KlbTLANO  Af  H AST  ) 


1 

wesr  0 HO I NANCE  Hi 

2 

*cI»UANO  ACCESS  Pf 

3 

PERIMETER  OH  SOOTH 

* 

perimeter  oh  west 

5 

PERIMETER  OH  NORTH 

6 

perimetfp  oh  east 

7 

CONNER  AVE 

a 

RIOGECREST  OH 

9 

H IHSCm  Oh  NORTH 

10 

AnThIS  AVE 

11 

OHAOShAM  ave 

12 

ELLIS  AVt 

1 J 

SAN  PABLO  ST 

U 

OAHLING  AVE 

15 

mALKER  ave 

lb 

f AIHCHUO  AVE 

17 

GENRIS  AVE 

18 

HIHSCH  OH  EAST 

19 

NEST  SANOIA  OR  TNI 

20 

PENNSYLVANIA  AVE 

21 

5JH0  ST 

22 

52NO  ST 

2 J 

50 T H st 

2b 

51ST  LOOP 

25 

h5TH  ST 

2o 

-6TM  LOOP 

27 

H9TH  LOOP 

28 

** 7T  H ST 

29 

-6TM  ST 

JO 

A ST 

Jt 

3** T H ST 

32 

NEST  SANOIA  OR  (*| 

33 

bOTH  ST 

3b 

30 Tm  ST 

35 

31ST  Pt 

36 

33R0  PL 

37 

35TH  PL 

38 

39TH  PL 

39 

-1ST  PL 

»0 

-3R0  PL 

bt 

3 7T  H PL 

<*2 

J2N0  PL 

»3 

38TH  PL 

bb 

*2N0  PL 

45 

■•bin  PL 

kb 
b7 
• 8 
*9 
50 
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»<i4N(  te  rs  t 04  nap 


KIHTIANO  AFb  4E  AS  1 » Hfn  MEXICO 


OH  IT  OP  >IK»FT  Li.  NOTH  Hi  A SURE  H£  Hi  = 0.00  INCH 

HAP  SC  AL  € * *00.  fHT  Pi  W INCH 

COON  J I n A T c SCALES  2.5  PLOTTER  I nChE  * PER  UNIT  * .189%  NILES  PER  COORDINATE  UNIT 

LtHuTn  CQH4  f RSI  OH  - .01000  HUES  PER  LENGTH  UNIT 

COOROlNATc.  USE  NOUE  *L  AST 

SCALE  RATIO  (TO  FIRST  HAP)  * 1.00000 

JcEAULT  SPEfJs  15.0  HPh 

OEEAULT  RATS*  2 

DEFAULT  SI3F>*  2 

DEFAULT  REFjSc  QUANTITY  ADJUSTMENT  FACTOR*  1.00 


l  10  4b  Q.Q  20/30  <9..1>  (9. .95) 

— - «*AO  JISTANCE  SPECIFICATION  ON  P PEVIOUS  LINE. 

The  HAP  COORulHATE  SCALE  ( 1.129)  OEV I AT£S  FROH  THE  OEFAULT  VALUE  ( 


HAPPLT  PARAMETERS  FOR  HAP  1 
SCR*  1.00000 
4NIN*  0.00000 
XL*  0.00000 
XSC*  2. >0000 
CNVLCN*  .01000 


T X*  0.00000 
XHAX*  10.00000 
XR=  10.00000 
TSC=  2.50000 


TT* 
THIN* 
VB  = 
P HOT* 


0.00000 

0.00000 

0.00000 

10.00000 


• 109  H1LES  PER  NAP  COCIROINATE  UNIT) 


YMAX*  12.00000 
VT*  12.00000 
PLEN*  70.00000 


2,2  410  5 2,2  420  5 ♦ 


2  20  27  0.0  JO  32  0.0  40 /JO  (.)  <5.85, .45) 

2 *0  32  0.0  >0/30  (.)  L 1 3 <5. .45) 

3 190  23  16.14  60  9 5,5  70  3 1,2  60  (6.3.1.15)  (4.5,1.15) 

4 AO  12  7.9  90  l 3 7.10  100  5 2,3  110  12  7.10  120  5 2.4  130/  (.1  LJl  (4,3.55) 

5 130  15  9.11  1*0  15  9.10  150  19  11. IS  160/<*>  (6.6,3.55) 

6 160  19  13.15  l70/(.)  (6.6.2.51 

6 170  22  I*. 17  160  / (.)  LS  (6. 3,1.4) 

6 1 A A 5 L.4  1 >u 

9 121  5 2.3  270  5 2.4  280  5 2.3  200  5 2.2  290  5 2t4  300  5 2,3  250  5 2,3  310  5 ♦ 

2,3  320  5 1.  « 330  / (. ) (6.35,3.3) 

1J  100  7 *.4  JbO  6 2,2  210  5 2.2  400  5 2,2  390  5 3,3  240  5 
2.2  *30  5 2.2  l 70  / (,)  L7 

7 140  5 1.2  200  15  10.9  210  11  8.7  450  3 2.0  470  7 6.4  220  5 3,2  70 
7 *50  6 0.5  *70  / (.)  LC 

A 50  * 1.1  r>0  5 2,1  230  22.5  15,16  2 40  15  9*10  250  5 2,2  150  4 1,1  260/  4.)  R26  ♦ 
(5.6.3.75) 

10  110  15  11.9  270  / (.)  R4 

1L  90  15  11.9  360  lb  10.10  260  / <•)  R4 

12  22 0 o 3.2  370  9 7,7  3A0  / (.)  Rb  (5.15,1.95) 

12  380  11  8,6  390  15  9.10  300/(«)  RA 

— - BAD  JISTANCE  SPECIFICATION  IN  PREVIOUS  LINE*  

HAP  DISTANCE  ( .111)  EXCEEOS  TOTAL  SEGMENT  LENGTH  ( .110). 

( 5.15  l.*S>  TO  ( 5.28  2.52) 


14 

290 

15  9.9  *00  12 

8.6  38  0 

/ 

(•)  R24 

15 

400 

75,3  >00  5 3,Z 

610  * 3,2 

180/  (5.5.1.45)  R2 

ib 

370 

5 2.2  230  6 3.2 

460  / ( 

• ) 

RC 

16 

480 

21  15.15  410  15 

9.9  310/1, 

) PI  9 

17 

500 

22  15.14  420  / 

(•)  LS 

17 

420 

15  9.9  320 

1A 

510 

22  l * . 1 5 430  / 

(.)  LS 

1A 

*30 

1>  10.10  330 

13 

100 

22  0.0  520 / ( . ) 

LS  (3*2. 

6) 

54 

520 

22  0.0  >30/35  ( 

.)  (3.3. 

75) 

20 

1170 

1 3A  0.0  2 60  31 

0.0  1240 

20 

0.0  530/25  (7.6  .3.7 5) 

20 

>30 

5 0.0  >*0  * 0.0 

545  5 0 

.0 

580  6 0,0  590  10  0.0  600  9 4.0  blO/25 

00000 
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***  «« Ss^^„PRwnciii, 


n>DD,C 
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GLOSSARY 


Air  Force  Refuse-Collection  Scheduling  Program:  a set  of  four  computer 
programs  that  perform  residential  refuse-collection  scheduling 
and  produce  printed  schedules  and  maps  of  the  routes. 

binary  search:  a procedure  for  finding  one  item  in  an  ordered  group  by 
repeatedly  halving  the  portion  of  the  group  that  contains  the 
item. 

free  format:  the  absence  of  card  column  restrictions  on  data  cards. 

map  coordinate  unit  (MCU):  the  length,  in  inches,  between  integral  divisions 
on  the  coordinate  system  appended  to  a map. 

map-description  data:  computer  card  input  that  describes  a street  map  to 
a computer. 

map  processing:  generating  numerical  data  that  enable  the  RCSP  to  produce 
an  approximate  copy  of  a street  map. 

node:  a numbered  point  on  a street  at  which  some  characteristic  of  the 
street  changes. 

parity:  a desired  modulo  2 sum  created  by  the  addition  of  a 1 or  0 bit  to 
a fixed-length  group  of  bits  during  transfer  of  data  from  core 
storage  to  tape  or  disk. 

pointer:  a variable  that  gives  the  location  of  some  other  variable. 

segment:  a portion  of  a street  between  two  nodes. 

shape  code:  characters,  either  two  letters  or  a letter  followed  by  a number, 
that  indicate  the  shape  of  a street  segment  or  string. 
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GLOSSARY  (Cone I'd.) 

spatial  clustering  of  streets:  selecting  streets  to  be  traversed  by  a 
► vehicle  on  one  trip  in  such  a way  that  the  streets  are  connected 

J by  other  streets  which  must  be  traversed. 

string:  one  or  more  connected  street  segments  having  the  same  shape, 
street  number,  speed  limit,  number  of  ways  of  travel,  and 
number  of  sides  serviced  on  one  pass. 
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